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Ya. Z. Tsypkin 


LENIN PRIZE LAUREATE 


On April 22, 1960, the Committee for the Lenin prize in the domain of science and technology for the 
Soviet Ministers of the Union of Soviet Socialist Republics awarded the Lenin prize to Doctor of Mechanics, Pro- 
fessor Yakov Zalmanovich Tsypkin for his work on the theory of sampled-data and relay automatic systems, pre - 
sented in the monographs, Transient and Steady-State Processes in Sampled-Data Circuits, Theory of Relay Sys- 
tems of Automatic Control, Theory of Sampled-Data Systems, published during the years 1951-1958, 


The editorial board of the journal Avtomatika i Telemekhanika (Automation and Remote Control) warm- 
ly congratulates Yakov Zalmanovich Tsypkin for his high honor, and wishes him further creative successes, 
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The V. I. Lenin prize for works on the theory of re- 
lay and sampled-data systems was awarded by the direc - 
tor of the laboratory of the Institute for Automation and 
Remote Control of the Academy of Sciences of the USSR 
to Professor Yakov Zalmanovich Tsypkin, member of the 
editorial board of the journal Avtomatika { Telemekhanika 
(Automation and Remote Control) and doctor of engineer- 
ing science, 

For the first time, a cycle of theoretical works in the 
domain of automation has received the highest award. 
The Institute for Automation and Remote Control of the 
AN SSSR (Academy of Sciences of the USSR) is rightfully 
proud that the works singled out for the Lenin prize were 
executed by one of its leading scientific workers, 

The editorial board of the journal Avtomatika i 
Telemekhanika (Automation and Remote Control) {s parti - 
cularly pleased to state that works of Ya, Z. Tsypkin, 
awarded this high honor, were published in this journal, 
starting in 1948, 

The bidgraphy of Lenin prize laureate Ya. Z. Tsypkin 
is a typical biography of a young Soviet scientist and 
patriot of our fatherland, 

Ya. Z. Tsypkin's first steps as a scientist were con- 
nected with the investigation of continuous control sys- 
tems, The concept of degree of stability, suggested by 
him in conjunction with P. V. Bromberg [1, 4, 5, 8, 18, 
48], entered into automatic control theory and was wide - 
ly disseminated in the international literature. His fol- 
lowing works were devoted to the dynamics of automatic 
control systems when lags were taken into account. Ya.Z. 
Tsypkin developed, and applied to the analysis of stabil- 
ity of these systems, frequency methods [4, 6, 7, 8]. The 
results which he obtained in this area underlay his doctoral 
dissertation, "Systems with lagging feedback,” [9]. 

The methods developed by. Ya, Z. Tsypkin for in- 
vestigating systems with lags received broad recognition, 
both in the Soviet Union and abroad, It is appropriate to 
recall here that, in evaluating the work of this period, 
Academician A, A, Andronov characterized Ya. Z. Tysp- 
kin as one of the most vigorous and promising scientists 
working in the domain of control theory in the USSR, 

The tremendous development of computing technolo- 
gy and radioelectronics in the postwar years engendered 
a heightened interest in discrete devices for control and 
communications, in particular, in relay, sampled-data 
and digital devices. Ya, Z, Tsypkin is one of the first 
scientists to have properly appreciated the importance of 
this trend and, more than ten years ago, he applied him- 
self to the development of a general theory of discrete 
automatic systems, His first work in this direction dates 
back to 1948 [10]. 

The fundamentals of sampled-data system theory 
were first presented by him in [12, 13, 15], The general- 
ity of the developed theories of sampled-data systems was 
not inferior to the generality of theories of continuous sys- 
tems, In 1951, Ya, Z., Tsypkin published a monograph on 
sampled-data circuits, Transient and Steady-State Pro- 
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cesses in Sampled-Data Circuit [17]. In [20] there was 
developed a frequency method of analyzing discontinuous 
control systems which allowed one to establish in which 
cases sampled-data control improves the system's dynamic 
properties as compared to continuous control, A further 
series of works were devoted to the effect of random stimu- 
li [22], of pulse shapes [33], of forcing actions [40] and of 
nonlinearities [23, 38] on the dynamic properties of sam- 
pled-data systems, 

Starting in 1956, he extended sampled-data system 
theory to cover systems with digital computing devices 
[39, 44, 46, 47, 49, 50, 51). 

The enumerated works were widely used by various 
organizations, designers, developers, investigators and 
users of sampled-data systems, The developed theory of 
sampled-data systems entered into the VUZ (university) 
curriculum and into courses for advanced engineering [45], 
It was widely used and further developed in a number of 
works by Soviet and foreign scientists, 

In 1958, Ya. Z. Tsypkin published the fundamental 
monograph, Theory of Sampled-Data Systems, [52] which 
is devoted to the analysis and synthesis of sampled-data 
systems and contains the fundamental general theory of 
sampled-data automatic systems, applied in various do- 
mains of technology. Together with sampled-data sys- 
tems, relay systems were investigated by him in a num- 
ber of works, In the theory of relay systems, Ya, Z. Tsyp- 
kin developed general and precise methods of investiga- 
tion, superficially similar to linear methods but, at the 
same time, applicable to this class of essentially nonlinear 
systems, These methods permitted the investigation of 
various modes of operation in relay systems, the comparison 
of types of circuits and the development of new, advanced 
schemes, These works [16, 19, 23, 27, 29, 31] were sum- 
marized in the monograph, Theory of Relay Systems of 
Automatic Control [34]. 

The monographs of Ya, Z. Tsypkin, and many of this 
papers, were published in the German Democratic Repub- 
lic, China, Poland, Rumania, Czechoslovakia, England, 
France, the German Federal Republic, the United States 
of America, and Japan, and were given as papers at inter- 
national congresses [43, 46, 50, 51). 

It is necessary to mention the characteristic traits of 
the creative style of Ya. Z. Tsypkin — clarity and explicit- 
ness in posing the problem, simplicity and elegance of the 
method of investigation, tense and concentrated purposeful- 
ness, 

It is no coincidence that his last monographs cons- 
titute three important successive steps in the solution of 
a single problem, that of creating a general theory of 
discrete automatic systems. 

Ya. Z. Tsypkin devoted much attention and work to 
the question of preparing a professional cadre, He is wide- 
ly known as a superb lecturer, His lectures on automatic 
control theory have been very popular [14, 21, 28, 45, 54). 

Under Ya. Z. Tsypkin‘s editorship, the fundamental 
foreign works on automatic control have been translated 
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and published here (McColl, Fundamental Theory of 
Servomechanism (1947); Lauer, Lesnick, and Matson, 
Servomechanism Fundamentals (1948); James, Nichols, 
and Phillips, Theory of Servomechanisms (1949); Flugge- 
Lotz, Phase Plane Method in Relay System Theory, 1959; 
and others). 

Under Ya, Z, Tsypkin’s direction, a large number of 
dissertations were executed, comprising a broad circle of 
questions in radio engineering, electronics and automatic 
control, 

Professor Ya, Z. Tsypkin is an important Soviet sci- 
entist, creator of a unified theory of sampled-data and 
relay automatic systems which has found wide practical 
application and has obtained world-wide recognition. 
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ANALYTIC CONTROLLER DESIGN, I 


A, M. Letov 
(Moscow) 


Translated from Avtomatika i Telemekhanika Vol, 21, PP. 561-568, May, 1960 


Original article submitted December 23, 1959 


The solution is presented for the problem of analytic controller design in correspondence with a given 
optimizing functional, The solution is given for a closed region which contains the boundary of excursion of 


the controlling organ. 


1. Optimal Systems, Defined in a Closed 
Region 

Essential in all automatic control systems is the lim- 
itation on the excursion of the controlling organs € : 





JEI<E (1.1) 


on the limitation on the magnitude of the velocity of this 
excursion 


1EI<f, (1.2) 


where & , f are positive numbers, 

With such limitations, the system becomes essential- 
ly nonlinear, since the presence of the equality sign in 
(1,1) and (1,2) permits the system to remain on the bound- 
ary of the region N of its definition, This region is closed, 
Moreover, it is impossible here to seek an optimal solu- 
tion to the problem among functions of the class C,, since 
there can be discontinuities on the boundary of the func- 
tion which solves the problem, By taking into account 
what has just been said, we formulate the following varia- 
tional problem. 

We consider a closed system to be controlled, in which 
the disturbed motion of the controlled object is given by 


th 
e equations (1.3) 


bx = te — (Dd) beans + mt) = 0 (k= i, ..., A). 


The notation in(1,3) has the same meaning as in [1], 
It is assumed that (1,3), in conjunction with the control- 
ler equations being sought, are defined in the closed re- 
gion N, characterized by (1.1), In this region there are 
given the sole boundary conditions 


Tro = H1 (0), «-~ 1 Mn(O) = Yno, &(0) = Fo; 
i (090) =...= tn(0o) =E(0o) = 0 (4) 


They mean simply this: that no matter what transient 
response arises in N, it must terminate, for t* = oo, with 
the system at the origin of coordinates, 

As the criterion of optimality, we take the integral 

co 
=| var, 


0 


(1.5) 


where V is the positively defined quadratic form 
V= Dia fp c&?. (1,6) 


We shall search for such continuous functions €, 14, 
. + «+ 2p in Class C (permitting, generally speaking, dis- 
continuities of the first derivatives) which minimize the 
integral in (1.5). 

The problem just formulated belongs to the class of 
the so-called discontinuous problems of the calculus of 
variations, To solve such problems there exist various 
methods of the classical variational calculus, as well as 
newer methods, presented in the works cited in [1], Here 
we shall use the methods of the classical variational cal- 
culus, supplemented by a nonlinear transformation which, 
in optimal problems of the first class [1], was usedin . 
[2}.° 
2. Solution of the Problem 





We consider the case when the limitation {s as in 
(1,1), and we set 
c= (0). (2.1) 
We define the function ¢ (€) as 


+& for t>°, 
9 ={e@ for [C1 <0, (2,2) 
—§ for C< —C, 


where ¢ * is a given positive number and ¢ (¢) is any 
continuous function whose derivative is continuous for 
|S} < &* , and which is equal to zero for€ = 0, In 
particular, such a condition is satisfied by the function 
y(&)= E sin € for ¢* =" /2 (Fig, 1). 

Transformations (2,1) and (2,2) translate the clased 
region N( E, 14,.. «» 1p) into the open region No ny, 
+ + «+ 1,) and permit one to use the well-known methods 
of variational calculus for the solution of ordinary Lagrange 
problems, The only special feature of this solution is that 
one must verify that the Weierstrass-Erdmann conditions 


* The author's attention was called to the transformation 
by I. A. Litovchenko, to whom the author expresses his 
gratitude, 
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hold at all points of discontinuity of the derivatives, To 
begin the solution of the problem, we set 


H=V + Dulin — Dorata ~ mp) - 
k 

















“ (2,3) 
eC) 
a 
-1/2 
x/2 C 
Fig. 1. 
On the basis of (2.3), we find that 
0H OH ’ 
ay a ao te Diets — >: ABes, 
on, k on, kVk — ak 
dH dH ff. , ] de 
=0 3 =[2020)— Drehe] FE - 
The equations of the problem have the form 
° - 
Te = Dy Yeats + muy (6), 
” a! 5 
hy =e b Kh +> 2a Gk 5 
dabahe + Bay (2.4) 


do 


Vs 


0O= [209 (c) — bY Maha 


They coincide to within the factor 9¢/9€ with the 
equations [(2.1) and (3,3)] of the analogous problem for 
the open region N(&, 1 ,... 4 pn) which was formulated 
in [1]. 

Therefore, in addition to the solution obtained in [1] 
for the open region, we should consider the new solution 
which corresponds to the equation 


. 5 


bal 
ue 








(2.5) 
By virtue of the definition given by (2.2), we set 
o(Q=+E fore (Sj >C. (2.6) 


This latter means that the optimal solution passes 
over the boundary of the region N (€, 14,..+5%p) 
For interior points of the region, the solution found has 


the form [1] , 
E= PaYa- 
Apa (2.7) 


Thus, the equation of the controller which corresponds 
to the functional of (1.5) is written as 





| —. 
> Pata for > Patjx | <5 
a a | 
”. a 2.8 
en g for >) Peta > 5, 
a 
"h7, E for >} Pata q —€. 





a 
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Obviously, this equation corresponds to an ideal non- 
linear controller with an infinitely large servomotor speed, 


The question of the stability of system (1,3), sup- 
plied with the controller defined by (2,8), should be con- 
sidered separately, It is the subject-matter of an inde- 
pendent problem. As was shown in [1], for a deviation 
|Zpaya | <§ , stability is guaranteed; for a deviation 
| Upa%a | >é a special investigation is required, lead- 
ing to the construction of the region of attraction, the 
boundary of which necessarily lies beyond the limits of 


| px %x|=§- Various cases can occur here, Thus, in 
example 1 of paper [1], where n = 1, by; = b, m, = m and, 
n,=1, we had the system of equations 

















4 = by + ms, 
» : kt+b _¢ 
| a gli tor Ps nm yj - ss 
“soa F |! nn 4] <5, (2.9) 
= for k+6 = 
=~ § = <-—é 
where the number k is defined by the formula 
fo , mea 
whale dhe aor a (2,10) 
é 
po (29) 
E 
0 ! 
j 7 
-E 
Fig. 2. 
é i E , 
—_—e—— ( j 





Sn 











Fig, 3. 


Figures 2 and 3 show the system's phase plane for the 
cases m, b < 0 and m < 0,b > 0, In the first case, the 
curve & =& (7) (2,9) is the sole integral curve by which 
the representative point returns to equilibrium for any 
No; in the second case, the system's stability is guaranteed 


mé 


Z* The features are analogous 


m> 0, b= 0. 


only for (%o)< — 


in the two other cases, in which 





n- 
2ed, 


and, 


2.9) 


10) 


- the 
ne 
nich 


inteed 


ous 





It is of interest to note that the choice of sufficiently 
small weight constants in the function of (1.5) permits 
an arbitrarily close approach to a relay characteristic of 
the optimal controller : 


' = Esign (—*+* 4). (2.11) 


However, the execution of the limiting transfer itself 
requires additional discussion, 


3. Possible Generalizations of the Func- 
tional of (1.5) 








We now consider the case of optimizing a functional 
of the form 


1G) = \ (V + &)dt. (3,1) 


0 


In accordance with transformation (2.1), we have 


t)'] dt. 





C Ta 1 tie 
I =(¢ @)=\ [Deenk + 09%) + (3 
a k 


The boundary conditions remain as before, i.e., 


(1.4). We set 
0¢ ¢)" 
og” 


H = Layyz + cg? (0) + ( 3 
+ a he |" - > OnxTa — Muy (| ' 
k a 











(3.2) 
Since 
0H 69\2; OH # so) 99 _ 
ae = 2(=) +a = [2c 70) +2 ahs a 


%, O¢ 
ou Dax ay * ° 
0% 
k 


then the equations of the variational problem will have 
the form 


a ' ie - 
. > Oka Ma + mg (8), 
a 


° 7 
y= — > baka + 2K Hx» (3,3) 


a 


se (B)']-f +2 








d 


Consider the last of these equations, After some ob- 
vious simplifications, we find that 


‘ do -'. oO? wo do 7 ayer 00 
2 [<2 te x : ‘| _— [2c¢ (0) — Yi mir | Fa (3,4) 
> » k > 


Therefore, (3.3) should be considered for the two 
Cases 


=o — Dra] o¢ : 


0 
K=O (3.5) 


when the system*s motion occurs on the boundary of re- 


gion N, and 
2 [SEE + SEC] = 2096) — mir 6) 


when the system's motion takes place on the interior 
points of region N, In the first case we have the solution 


e(S)=té ior» [SI >0 (3,7) 
In the second case, since _o9 F on # 9(f) E q.(3.6) 
=, ’ So ee 


reduces to the form 
- P a ‘ 
29 = 2cm — >) mex. (3.8) 
k 


In conjunction with the first 2n equations of (3,3), 
this determines the controller equation 


F = >) Paha — re. (3.9) 


It is essential to note that, inasmuch as E qs.(3.3) of the 
variational problem coincide exactly with (2,1) and (3,3) 
of [1], the equation just found, (3,9), coincides exactly 
with (4,5) of [1). 


On the £,4,..., np phase space we sketch the 
two hyperplanes 


2) Pate = 7%. (3,10) 


Then, the controller equation being sought will be 
written as 


> - =k 


+- E tar >: .? 





| am 


} r 





: = LPaNa _ re for 


for 2Pa%q 
r 








svrl 


4, The Weierstrass-Erdmann Conditions 





To complete the investigation of the necessary con- 
ditions for the existence of an optima! solution of the 
analytic design problem, one should convince oneself 
that the Weierstrass-Erdmann conditions hold at the nodal 
points of function €. 

In the case when one is considering a functional of 
the form of (3,1), these conditions reduce to the holding 
of the following equalities; 
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0H*+  dH- dH+ _ 0H- (4.1) 
mu, mM & (hm ft, .. + Be 
a A P ‘@H +-\>- / OH 
H* — (Se "—(4t) =a p( aes) (St), a. 
~ an, a / 2 (a3, 7) a 


Here, H is the function in (2,3) or in(3,2), and the 
"+" and "-® signs denote that the computation of the 
corresponding functions is carried out to the left or to 
the right of the nodal points of function &. 
One easily convinces himself that, if the first n con- 
ditions of (4,1) hold, which reduce to the equations 
at = ie 


(k=1, 00-5 2)» (4,3) 


then the two remaining conditions are also satisfied for 
any functions §,4,...,5 Nn of class C, As for (4,3), 
it can always be satisfied, as is easily verified, if the op- 
timal system is stable, 

Thus, in the example considered in Section 2, the 
solutions for the optimal controller in the case when 
have the form |¢|>& 


n= c,e% 
= ce-¥ a 4 ae oa. ott, (4.4) 
and, in the case when |€|<€ 
eee 
~ coe [— (k +b) cye-*t + (kt — b) Cye**], fit 


pa NH e+ ee + (bb) eet. 


From the conditions on the left end (for t = 0) 


— m < = 
to = FE B= ts. (4.6) 


From the conditions on the right end (for t = o), 
Co = 0. (4.7) 
The continuity conditions at the nodal points are 


- m k+obd- 
Cy = cya —_ = c= +6 





(4.8) 
The single condition of (4,3) has the form 


2am— acy 


— bt, 
Coe fs he -}- —_ 





ete = — = (k + b) cy. 
(4,9) 
The last equation of (4,1) and (4,2) reduce to ident- 
ities, 
Equations (4,6), (4.7), and (4.9) define the constants 
Cy, C7, aNd Cy; (4,8) define c, and the time te, corre- 
sponding to the nodal point of function §. We, hence, 
find that 
tL 


ehitinns ni 
- =(%o+ =| ete, (4,10) 


om 
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It is easily established that, form < 0 and b < 0, 
(4,10) is solvable for t® >0 for any Ino | > | . 
if this latter inequality does not hold, the system does 
not go beyond the boundary /|€| = € and no nodal 
point exists, In the case when m< 0 andb > 0, the 
equation has a solution with respect to t* > 0 for any Nq 
which satisfy the inequality 


~S 


- _ < to|S 


k+b 





This is againthecase shown on Figs, 2 and 3, The 
analogous features also uccur for the other combina - 
tions of signs of m and b, 

If the functional of | (3.1) is used as the criterion of 
system optimality, then the Weierstrass-Erdmann condi- 
tions also are met, Indeed, the solution of the system for 





|E| = is as follows: 
<= mE - 
n=aetee, — = §, 
(4,11) 
h = ce it - ee. 


With the boundary conditions at infinity taken into 
account, the solution of the system for || <t is 


4 = ce"! + cyeHst, 


e_ m—95- wt 4. ta om 
7 od m cael, (4,12) 


ee 
ve+6 -* 


We find from the initial conditions that 


— mE ~ = 
jYo= at, &=€. 


7 (4,13) 


The continuity conditions for the functions § and 7 
at the nodal points are 
ie b > oo b s 


= ———e, + C. 
m ‘— m - 


— = <S 


bt, - ! 
Cee FH 1 + Ce. 


(4,14) 
The Weierstrass-Erdmann conditions reduce to the 
holding of the two equalities: 


von, (EY = (ey 


or, what amounts to the same thing, to the equalities: 


evry 





(4,15) 


E - ~ 
mas 2amE 4 ae pdt, — 2a 2a 


oe Sse aa eee 


coe! 


a (t41 — 6) Cr + pe (Ye — b) Co = 0. 





or 


3) 





Equation (4,13) defines the constant cys from (4,14) 
we find that 


@ 1 ‘ - 
Ee ae (42 — 5) n° — mf}, 


(4,17) 





(vE — (ur — 5) 1], 


ae 
a ta Ha 


where n* =n(t.). The first equation of (4,16) defines the 
constant cy, the second defines the moment of time te 
when the system passes through the nodal point, This 


last equation, in correspondence with (4,17), takes the 
form 


mE (py + pg — 6) = (1 — 5) (pe — 5) ie (4,18) 

This last equation means that, at the nodal point, the 
servomotor speed must vanish, as was also the case with 
(5.4) of [1]. 
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The problem of finding the optimal control law in pulse-relay second-order systems is considered. 


Let there. be given a physical system which is de- 
scribed by the n-dimensional vector x;(t) and let the posi- 
tion of the representative point in n-dimensional phase 
space at the initial moment of time be defined by x;(0). 
The representative point must be translated to the given 
point €,. The process will have been completed after 
time T, and will be optimal if no other process provides 
for the attainment of €; in a time less than T, 

The optimal control problem for systems described 
by differential equations of the form 
dz; Q) 

ar = Pi (Ty Ur, L) UY anO 8... +, M be 8, 1,2--6 ¥) 


has been widely investigated in a number of works [1, 2]. 

In this posing of the problem, time is continuous, and 
one imposes limitations of the following type on the con- 
trolling quantities u r 


(2) 


In [3] the analogous problem .is solved for systems 
described by finite difference equations 


}uy|< Ni. 


x; (ke + A] kh) — 24 (kh) = [> asx; (kh) + 
j=1 


m 
* Zz git; (kh) +-f; (kh) |h. (3) 
l=1 
In this problem, time is discrete, t = kh (k = 0,1,...) 
for an arbitrarily chosen h, and conditions of the follow - 
ing type are imposed on the controlling quantities: 


| w (kh) | <N. 


Such a system was investigated in [4]. 

The system of second-order finite difference equa- 
tions which we consider in this paper, these equations de - 
scribing planar motion, has the form 


x({k + 1) h) — x(kh) = y (kh) h, 
y ((k + 1) hk) — y (kh) = u(kh)yh, 
where x, y are the point coordinates on the phase plane 


and u is the point acceleration (controlling action), The 
behavior of the points is considered under the limitation 


(4) 


(5) 


(6) 


| ux | = const, 
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In this case (a pulse-relay system), it is required to 
determine the control law { m4 } which allows the initial 
point (X9, yg) to be translated to the origin of coordinates 
of phase space in the least possible number of steps K®, 
The time here, as in (3), is discrete but, in contradistinc - 
tion to (1) and to (3), the values of the controlling quant- 
ities are also discrete. 

The optimal control problem for a system with (5), 
which we shall solve in the present work, can be classi- 
fied as a simple number-theoretical problem. It should 
be mentioned that problems of this type will continue to 
arise with ever greater frequency in connection with the 
use of discrete computer technology in automatic control, 

It seems advisable to sketch briefly our plan for find- 
ing the optimal control law for such a problem. We first 
find the optimal control law for initia! points of the form 
(n’, 0), for which we shall use the results of the work of 
N. N. Krasovskii [3] (the corresponding expressions for the 
number of steps K® and for the control law { u’} can be 
easily and immediately obtained). We then construct the 
optimal control for all the integral points of the x axis, 

x > 0, and we show that the optimal trajectories completely 
fill up the region of the (x,y) plane defined by the expressions 
0 =x = y(y-1)/2, y< 0. Finally, we shall show that the 
optimal trajectories of all the other integral points of the 
(x, y) plane pass through one of the points already considered 
on the x axis in one of the half planes x > 0 or x < 0. On 
the basis of the optimal control law obtained, we shall 
construct an optimal switching line. At the end of this 
paper we shall discuss questions related to the nonsingle- 
valuedness of the control law obtained. 

We turn now to the solution of the problem, to be 
carried out in accordance with the plan just given, 

The solution obtained in [3] for a second-order sys- 
tem with (4) renders it possible, for a given h, to choose 
an optimal value of | u, | which depends on xj» and &j, 
but does not depend on the number k of steps, and to de~ 
termine the unique moment for changing the sign of the 
controlling action (corresponding to a transition from a 
type N to a P-trajectory [1]). 

The peculiar feature of our problem, when (4) been 
replaced by (6), is that, generally speaking, for a con- 
stant modulus of point acceleration, the sign of this ac- 
celeration changes more than once during the control 
process, 








Indeed, we consider a point lying on the x axis ata 
distance X» from the origin of coordinates. Let x and y 
be changed for integer steps, i.e., | Uy, | si,h=1, if 
|u| =a and h =b, where a, b* 1, this leads only to a 
corresponding change in the scales of the process, In ac- 
cordance with [3], when (4) holds, K° is defined as the 
least integer which satisfies the inequality 


K2h? 
| 2%. 


Replacement of (4) by (6) entails only the change of 
the inequality sign to an equality sign, i.e., 
K*h" 


“er = Io. 


(7) 

In order that (7) have integer solutions in K, it is 
necessary that xy be a perfect square, i.e., x)= n*. For 
such points of the x axis, the optimal control law of [3] re- 
mains valid, it taking, in the given case, the form 


u® (kh) = — sign [> —k| ’ 


K® = 2V x = 2n. (8) 

In the case of (6), for all the remaining points of the 
x axis, it is not possible to reach the origin of coordinates 
with a single switching of the acceleration's sign. 

We limit ourselves initially to a consideration, for 
such points, of the half-plane x > 0, and to the require- 
ment that one attain in the least number of steps some 
neighborhood of the origin of coordinates consisting, for 
example, of the two points (0,0) and (1,0). We define 
the position of the point on the x axis by the distance p: 


Ty = n* + p, 
where 


O<p<2n+1. 


Since the chosen neighborhood is characterized by 


.¥@ = 0 then, for yg = 0, we have the equation 


yu=— J yy 


uj >0 uj<o 


(i, 7 =1,2,..., K%, 


(9) 
i.e,, the neighborhood of zero may be reached only in 
an even number of steps, It can be shown that 


p=1 
p=2,3,...,2n+1. 
(10) 

In fact, for a zero initial velocity, the path taken by 
the point in 2n steps is maximal for the choice of ug(kh), 
according to (8), and is equal to n*, With this, the point 
is at the distance p from the origin of coordinates, so that, 
for p= 1, the point falls in the neighborhood of the origin 
of coordinates, and (8) defines its optimal control law. 


When (9) is taken into account, for p =* 1, the magnitude 
of K® cannot be less than 2n + 2, 


2n, if 
~ (Q2n+2, if 


"0 


We now assume that the optimal control law has the 
form: 


u, = — sign| 4 —k— 1] — 28. 
ji for z=k 


é.- = 
2 \0 for zk (k=4,2,3,..., K%, 
where K° is defined by (10), 


(11) 


Now we show that one can choose a z corresponding 
to anoptimal control.* As the end point for an even p = 
= 21, we choose the point x,0= 1; forp=21 +1, let 
such a point be xx = 0(2= 1, 2,...), Then z must be 
so chosen that, for any parity of p(p* 1), ’the following 
equation holds: 

se + B(f]—e [=e + 
where S; is the path taken by the system during the con- 
trol process in Kg steps; E(a«) denotes the integral part of 
a, At the same time, by adding the segments of path 
taken by the system for the controlling actions of (11),we 
obtain 


Sxe = n? + 2[2n — 1—(z — 1)] — 1. 


By comparing the two last formulas, we obtain the 
expression for z: 


z= K°—E|PS*). (12) 


Consequently, if z is defined by (12), then (11) defines 
the optimal control law. 

The optimal law just obtained is easily extended to 
any integral point of the phase half plane x > 0. In fact, 
the initial points x» = sign yg for x» = yo(y9-1)/ 2 (the 
boundary of this region is the N-trajectory which terminates 
in the origin of coordinates, {.e., X9 = yg(Vo-1)/ 2) lie on 
one of the optimal trajectories considered, One easily 
convinces himself that at least one optimal trajectory 
passes through each point of this region. This follows 
from the fact that the optimal control law, up until the 
first switchings (i.e., for ks K°/ 2-1), does not depend 
on p. The points sign xX» = sign yp and sign x9 = ~sign yo, 
with the condition that x»< yg(yg-1)/ 2, are taken over 
into the case just considered after y» steps u, = — sign X» 
in the right-and left-hand planes respectively. It is easily 
seen that the “initial translation” of points moving along 
P- or N-trajectories are also optimal. If necessary, the 
optimal control law can be written for arbitrary initial 
conditions, For example, for the case of initial points 
defined by the conditions 


sign % = —signyo, ZT>Veyo(yo—1) (p+ 1), 
K® = 2(n + 1) —| yo}, 
uk = — sign [n — | yo| —k] — 28,., 


* We adopt the convention that sign 0 = 1, 
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we have 


z= A — Ble —|Yol, 


P = % + */2Yo(Yo + 1) — n°. 


The control law just found corresponds to a transla- 
tion of points to the neighborhood of zero, but it is easily 
verified that if K° remains invariant, and for any number 
of switchings, if (9) is met then the parity of S; is un- 
changed, Therefore, the law that translates the point di- 
rectly to the origin of coordinates must be supplemented 
by two steps (when u,o, = - UXO, = 1) to all the tra- 
jectories which terminate in the point (1,0), Then, the 
optimal control law, in the form of (11 and (12), will be 
vaild for all points and, for all p, we shall have K®* = 2n+2, 

The control law in the form of (11) and (12) allows 
us to find the optimal switching line. 

We consider the phase half-plane x > 0 and the half- 
line y= y®* on it, It is obvious that the optimal stimu- 
lus for all points with integral coefficients on the half-line 
y=y°* > 0 will be u® = ~sign x = -1, Now let y® <0 
and v=| y4 (Fig. 1). The line under consideration inter- 
sects the N-trajectory x = v(v-1)/ 2at point { af wv-1)/2, 
y*}. Then, for all points with coordinates 0 =x = 
< v(v-1)/ 2; y = y®, we have that u’ = sign x= +1, It 











y 
0 y? (ust)? 
Ls a re 
a“ 
a ao -s 
Fig, 1. 


remains to consider the points x > v(v-1)/ 2, y = y*>0. 
We select an interval of the x axis which contains 2v + 1 
points v? < x, < (v + 1)*%. For all points, according to 
(10), K°= 2v7+2 (in the sense of a translation to the 
neighborhood of the origin of coordinates) and (11), it is 
required that u° = sign x = 1, since a point of such a seg- 
ment falls on the line after v steps, For points on the x 
axis such that xg = (v + 1)*, the first change of sign occurs 
on the step whose ordinal number is known to be greater 
than v, and for these points u® = ~sign x, Thus, for yg = -v, 
the optimal control is these points u° = - sign x. Thus, 
for yg = —v, the optimal control is 


r>b, 
r<b. 


wu? = 


—sign7y, if 
+-sign7,. if 


The abscissa of point b is obviously equal to 
v(v—1)/2+(2v+1)=(v+ 1)(v+ 2)/2. Since the argument did 
not depend on the choice of v, one may assert that the ex- 
pression just found defines the optimal switching line on 
the half-plane x > 0, 
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By carrying out an analogous argument for the half- 
plane x < 0, we can obtain the complete expression for 
the optimal switching line: 


waft if e>s(e 2), (13a) 
{- l, if r< f(r, zt), 
where v = ly| ® 
Jae) _ EERE TA. if 2z>0, 
I(v, 2) = v(v + 3) 
lf. (v) = —. — 9 if zx < 0. 
(13b) 


In the half-plane x < 0, to retain the "2" and * <* 
signs unchanged, the switching line must be translated by 
unity toward the y axis, 

The switching line just obtained is optimal in the 
sense that the representative point is translated either to 
the origin of coordinates or to a neighborhood of it. 

The optimal control law, which we have found in 
the form of (11) and (12), is not unique. We now show 
that it is possible to define other forms of the control law, 
We consider the control of points on an expanse of 2n+ 2 
steps, A control for which the sign of the stimulus is con- 
stant (for example, a negative stimulus) on the first n + 1 
steps (i= 1,2,...,m+1) and, fori=n+2,..., 2n+ 2, 
the sign is reversed, we shall call a “canonical” control, 
With canonical control, there are (n + 1)? points on the 
path, Since, according to (9), the numbers of steps with 
"+" and with "-" signs are equal, then all possible con- 
trols can be obtained from the canonical one by per- 
mutations of signs pairwise between both sign-invariant 
portions, If, for example, we change the signs of the two 
terms with ordinal numbers m, and m, 


1<m<n+1,n4+2<m,< 2n-+ 2, 


then the points passed through on the path will be fewer 
than with the canonical control, and will be defined by the 
expression 

Sento = (n -}- 1)? —2(m, 


— my) (14) 


If the initial point is defined by the conditions x = 
= n*4 Pp. y = 0, then we can obtain the expression for its 
"optimal" control laws by equating Syn +, =n? + p. From 
this, we obtain the relationship for m , and mg: 
2n—p 
iia (15) 





M,g— mM, = 


| | 5 16 17 10 
pt EUG 


Fig. 2. 

















In this work, we fixed the value of m, = n + 1, and 
my varied as a function of p, Such an approach made it 
possible to obtain the optimal switching line. 






































A typical form of controlling stimulus is illustrated 7 u* y 
on Fig. 2, Here, Xg= —19 and yy = 0, System's sig 
Figure 3 shows the optimal trajectory for the point R, , R, 
(Xp = 12, yo = 0) which corresponds to n = 3 and p = 3, 
Trajectory 1 is formed in accordance with optimal switch- Py ; 
ing line 3, and trajectory 2 (dashed lines) is obtained from if) 2 
the control of (15) for m, = n, Fig. 4, 
y 
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THE EFFECT OF RANDOM NOISE ON THE STEADY-STATE OPERATION 
OF A STEP-TYPE EXTREMAL SYSTEM FOR AN OBJECT WITH 
A PARABOLIC CHARACTERISTIC 


T. I, Tovstukha 
(Moscow) 


Translated from Avtomatika i Telemekhanika Vol. 21, No. 5, pp. 575-584, May, 1960 
Original article submitted October 21, 1959 


The optimal parameters of the controlling portion of the extremal system under consideration are calculated. 
These parameters must be determined from estimates of the deviations from the minimum characteristics in the 
steady state, 

It is shown that, in the particular case of uniformly distributed noise increments, it is most advantageous 
to operate with small steps, to be determined below. 















































1. Posing of the Problem were 7} 
ig ! . 2(n) 
In this paper we consider the case of a step-type ex- rin) | sy, yi 7 
, “+ 1 
tremal system whose object has a parabolic characteristic ! z 1 juin) 
y = x*, The system's task is the automatic search for, and an ------- , 
maintenance of, the minimum of the characteristic. The 
following posing of the problem, as well as the notation, < 
are the same as in [1] of A. A, Fel’dbaum, All quantities Fig. 1 


are considered at discrete moments of time t = nT (n= 0, 
1, 2,...), where T is the time between two successive 

working steps, or the duration of one cycle, Figure 1 is a 
block schematic of the extremal system under considera - 
tion; O is the object of control, its input and output being 


operation which implements this motion. The controlling 
portion, with its testing motions, disturbs the object's 

regimen. In the system under consideration, at the beginning 
of each half cycle, the object's input x [n] receives a step- 


related by the function wise increment of (Ax), = a (at the beginning of the cycle) 
and (Ax), = -a (in the middle of the cycle), At the object's 
y(n] = 2? [n}. (1) output, at the end of each half cycle, one measures the 
Random noise z[n] acts at the object's output. The corresponding disturbances w, = yy — and > = and 
noise increments v = z[n + 1] -z[n] are assumed to be compares them, Thus, one determines the direction of mo- 
statistically independent quantities with probability den- 40M toward a alae, or the sign of the working sep (the 
sity p(v) and zero mean value, The noise z[n] can be magnitude of the working step, as of the test step, is a): 


identified, for example, with the random noise in the Az[n] = —asign Aw[n], 


measuring block of an actual system. Slow shifts of the itd «eo (8) 
characteristic, or input noise of the object, will not be is “ ” — a 
considered here, The action of noise retards the automat- It is possible with this to make a . oo _ , a de 
ic search for the minimum of the object's characteristic, ing away from the minimum, since the quantity w Is 
where the minimum itself is determined with a certain up of the true output of the object plus noise, If one know 
exser the distribution law of the noise increments, p(v), then 

The actual object output, which can be measured, is  “€ probability of a false step from the point x equals 
the sum of the output y[n ] and the random noise, i.e., the p(z) = P{v> Ay}, 
quantity where Ay = yy~ Ya, V = 247Zg, OF 

w[n] = y(n] +2[n]. 9 tool 
vs piz)= | ple)de. ( 

The quantity w[n] is applied to the input of the con- —co 

trolling portion of the system ,C. The probability of a correct step (one directed 


The search for a minimum is implemented by means _ toward the minimum) from point x equals 
of successive cycles; each cycle includes a test to determine 
the direction of motion toward a minimum, and a working q(x) = 1— p(2). 
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Possible neighboring points x are at a distance equal 
to the working step a from each other, We renumber the 
possible discrete states of the system, assigning the sub- 
script zero to the state which corresponds to the minimum 
of the characteristic, i.e., 


r= ta oog—1,0,1,...,00), (5) 
We introduce the notation: Py = p(x;) fs the probabil- 
ity of a false step from xj, qj = q(x;) is the probability of 
a correct step from x;. 
If pj, {8 the transition probability from point x; to 


point x), then, in our case, 


(i = —oo,, 


={? (k>0), 
Peek+i = I (k <0), 
- (k > 0), 
me ls oe 
(6) 
Pri =O (lst k —1, k + 4). 


Thus, the process of seeking a minimum can be con- 
sidered as a discrete Markov process, 

We denote by the "weight" of state x, the mean value 
of the object's outputs for the two tests, using the notation 
(Ay), for this quantity, In correspondence with (1) and 
(5), 


(Ay: =a2(?+1) (i =O4442%-.- 7 


The difference of the values of y in the test, Ay = 
= y; = Yq, is defined as 


(Ay); = 4ia? (§=0,+4,42,...). (8) 


The probabilities of false and of correct steps from 
state x; are 
—4\i|a* 


r= | ple)de, 


—oo 


R= 1— pi. (9) 


The steps from the minimum, both at the point x, =a 
and at the point x_, = ~a are false, with a probability 
of pp = 1/2, 

For convenience and clarity of the exposition, we 
can construct an equivalent scheme of the discrete Markov 
process [2]. Such a scheme is shown in Fig. 2. The points 


State 
— 


; . weight 
fn 7 
a 


(1) (2) + State 
ordinal 
number 





Possible states 


Fig. 2 


on it represent the possible states of the system. ‘The ar- 
rows show the transitions from state to state, The'prob- 
abilities of these transitions are shown, Z is the operator 
for translation in time, showing that the transition’ is car- 
ried out during a working cycle, 


We denote by P,[n] the probability that the system is 
in state x; at time th = nT, Then, the mathematical ex- 
pectation of the output quantity is 


co co (10) 
Maly) > > (Ay) Pitn] =a® >) (@ + 4) Pafn). 
=—0O 1=—0o 


It is necessary to determine lim M,[y], if this limit 
Ir? 


exists, and also to find the a for which this limit is a 
minimum, 

The case of an extremal control system with the 
piecewise -linear characteristic y = | x| was devoted in 
work [1], It was shown there that, in the steady state, 
the minimum value of the mean error in y corresponds 
to a step of size a © 0 for normally or uniformly dis- 
tributed output noise, 


2. Determination of the Mathematical 
Expectation and the Dispersion of the Ob- 
ject*’s Output Variable in the Steady State 











We write the equations for the probabilities of the 
system's states; 


P_, [n}] = P_s[n — 1]q3 + P_, [n— 1] py, 
P_,(n] = P_(n—1]}q2-+ Poln—1)4, 


Po{n] = P_,(n— 1} + Prin — 1] 91, (11) 
P,|[n}] = Pa[n—1]q2 + Pol[n— 1] - 


P,(n] = P3[n— 1] qs 


+ P,(n— 14] py, 
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Fig. 3 


We have taken into account here that, in correspond- 
ence with (9), p; = pi; and q; = qi. 
We now consider the finite system of 2k + 1 states, 


shown on Fig, 3. The equations of such system are writ- 
ten in the form 


P_, [n] = Piss [n — 1) Pr—i + P_, [n — 1) Pk» 
P_x41[n] — Perse [nm — 1] Pro + P_e(n — 1] Ge: 


P_,[n] = P_,[n— 1] p; 

P_, [n] = Po[n — 1] > 
Po[n] = P_,[n— 1] q 
P,[n] = Py[n — 1] > 
P,[n] = P,[n— 1] p 


+ P_s[n— 1} Qa 
+ Pi[n— 1} q, 
+ P[n— 1] qs, 
+ Ps[n— 1}qs, 


.. 0 £4. © #6 #8 12 @ 9's 2 Oe 


(12) 
We consider the limiting system for n- oo, assuming 


that the lim P;[n] = ry exists, In this case, 
noo 


Toe = T—e4+i1Pr—1 + —KPr» 


Tk = T—k42Pa—2 + T—KIk> 
rg =T-1Pi + 7-393» 
r4=lo= Pde» 
1 > + 7-242 (13) 
To = 7-191 +7191, 
4 
n=loy + TeG2, 
Te =T;P1 +1393» 
Te—-1 = Th—-2Pr—2 + TAK» 
Tk =Te—1Pr—1 + TkPk- 
It is obvious that ‘ 
> y= | 
i=—k (14) 
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We always have 2k + 2 equations in 2k + 1 unknowns, 
We consider the system of the first 2k equations of (13) 
and (14), and we determine the rj which satisfy this sys- 
tem, It will easily be shown that they also satisfy the 
last equation of (13). 

If we also take into account that, by virtue of sym- 
metry, tj = r-j, we shall consequently consider the sys- 


tem 
Po = 2ri91; 


1 
M1 = Toy + 292 


Te = 71P; + TsQs» 


1 
Pe—1 = Te—aPr—2 + TeTk (15) 
k 
rot+2> n=1. 
i=1 


We shall successively express all the r, in terms of 


To = 29171; 
1 
=n 
Po 1 
bp. 1 gs 4s 


It is easily shown that, in general, 


Ti =TyPydghs. . . N—G (2<i<hk), 16) 
where 

he one Cn = 2 (m =1, 2 k) 17 

m Im Filial (17) 


By substituting, in the last equation of (15), the ex- 
pression for rg from the first equation and the expression 
for rj from (16), we obtain 


1 


k 
2 (1 a qi oa > Prideds see d_i8s) 
i= 2 


The value of ry is defined in terms ofr, from the first 
equation of (15) 


n= 





(18) 


To = 27191. (19) 





16) 


17) 


on 


18) 


(19) 





It is easily verified that the thus determined rj satisfy 
the last equation of (13), f.e., the one we didn't consider. 
Thus, the state probabilities in’the steady state are com- 
pletely given by (16)-(19). 


We now turn back to (10), We need to determine 
the magnitude of the mathematical expectation in the 
steady state, f.e., 


lim M,[y] =a? >) (i? + 1) 7. (20) 


i=—oo 


oo 


In (16), (18), and (19) we let k tend to infinity, and 
we substitute the values of r; thus obtained in (20): 


co 
1+ Cit >) (i+ A) Adra. «Ayah 
i=2 


= q? 





oo * (21) 
i +- Ci -{- >} Ayhs eee Ayi& 
i=2 


Knowing the rj, we can easily determine the disper- 
sion of the quantity y: 


Diy}= >) (Ayn — M* ly] = 


t=—0oo 


co 
Cit >) ifhida... Ay_it 
4 i=2 
=a 





co 
i++ bY AyAg.. > Van A 


i=2 


3. The Case of a Uniform Distribution 
Law for p(v) 








In this particular case, the computations of M[y] and 
D[y] by (21) and (22) are easily reduced to numerical re- 
sults, Let (cf., Fig. 4, a) 

1 ' 
«I? I<), 
(0 (jel >). 


In correspondence with (9), the probabilities of false 
and correct steps are found from the formulas 


p(v) = (23) 


1 F 4a* /. n 
=" (i< Gat}? 


R= 
. (i> qa) (24) 
sti (i< qi), 


a= (i> zs) (i =0,1,...). 


We denote by k the integral part of the expression 
(n~4a*)/4a®, n = 4a" + B, y = B/ 4a* (Fig, 4,b) 
Then, from (21), 


k+1 
1+ Cit >) (iP + 1) Ade... Af; 
M [y] = a* we ’ 
1+Cit >) Asds...A, 4% 
t =2 





where — and €_ are determined from (17), with ac- 
count being taken of (24), 














+2 Aly) _ 4) marty). (22) 
AY) 
a 
-9 0 7 v 
Flu) 





wart TTL 





JN ‘0 D 
¥ ~tne? -4a? 4 
Fig. 4 


To simplify the computations, we shall consider, 
not all a, but only those for which 
n 
las (25) 
where m is an integer, 


we easily obtain from (20), (23), and (24) (cf, the 
appendix) 


eta (Z>1), 


‘ni ) 
ae 1.5a? (<1). 7 


Here, since M[y] is a continuous function of a, n/ 4a” 


can take any value, either integral or fractional, 


We compute the dispersion Diy] analogously (cf., 
the appendix) 


: @eeoe 
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or we compute the mean-square deviation 


VaG-*) Gr) ag 
oly)= a (2<!) 
2 4a* */}* 


The dependence of M[y] and 9[y] on a for several 
different values of n is given on Fig. 5. For a comparison 
on the same graph, we give the dependence of the math- 
ematical expectation and the mean-square dispersion on 
the size of the working step for the case of the piecewise - 
linear characteristic y= |x|, considered in [1], for the 
same noise distribution law. The corresponding computa - 
tions are given in the appendix. 


Q 
0 





Fig. 5. Comparison of the curves (M/n) (a/n) 
(solid lines) and ( o/n) (a/n) (dashed lines) for 
parabolic characteristic and modulus character - 
istic in the case of a uniform distribution: curves 
1 and 5 are for a parabola with 4 = 4, curves 2 
and 6 are for a parabola with n = 1, curves 3 and 
7 are for a parabola with n = 0,25, curves 4 and 
8 are for the modulus. 


smallness of the mathematical expectation and the dis- 


of the step decreases the mean-square deviation insignif- 
icantly, but increases the mathematical expectation of 
the output, In the modulus case, y =| x| , it is better to 
make the step larger than for the relationship y = x” 
since, for the modulus, an increase in the step (for small 
values of the step a) gives an insignificant increase of 
M[y] accompanied by a decrease of the dispersion, 


In conclusion, the author wishes to thank A, A, 
Fel'dbaum for posing the given problem and for his use- 
ful advice, Thanks also to S, Ya. Raevskii for discussing 
the results obtained, 
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It is clear from the graphs that, for the characteristic 
y= x’, the optimal (from the point of view of simultaneous 


persion) working step a is close toa = 0, A small increase 


APPENDIX 


1, Determination of the Mathematical 
Expectation and the Dispersion of the 
Output Variable for Uniformly Dis- 
tributed Noise 














By substituting (25) in (24) we get 


1 i | “eae 
Se (i< m), = {tt im (i < m), 
(0 (ti >m); \1 (i > m) 
(D 


We find from (21), for m > 1, 


m 
1+ 21 -+ >) (i+ 4) Agds...Aj_.f; 
i=2 
M |y| =a? — , 
’ ~ ¢ 
1+ti+ » Aghg... Ay G; 


i=2 





where A; and © ,are determined from (17) with account 
being taken of (I). 
After some transformations, we obtain the expression 

















25) (i244) mi antl 
e. (m — i)!(m + i)! 
M [y] = 2 _ m 
ae 
S (m — Yim +! 
™m 
We let A ee . Then, sin 
(™)= 4 iH—jim+i!* eis 
1=0 
2m—1 { 





2m 

y - (2m)! 

= 72m — 7)! ~ 
j=0 


2™ then A(m) = lim)! + se 


“ 


and the denominator of M[y] is written as 


22™m '2 


_ (2m)! are 


Ee om 1 
amP = leat + =n st 


We now consider the numerator of M[y} 





es ’ m!2 
23+ Veep! = 





i=o 
— Im!2(m? a 
= 2m! (m } 1) Pa (m — iylqm ; ip! = 
i=o 
™m 
‘ m2 — i2 22M™—1 p12 
a 12 pee = Ps 8 Le 
a 4 mh ‘-“—_ 
i=0 
Whereupon M[y]= a%(m + 2)/2 form > 1, 


By taking into account that m=  /4a?, we obtain 
finally, 
> i. 


* ~ 
M |y|= = + a* for ia? > 









n), 


n). 


unt 


ssion 


since 


ybtain 








integer m less than unity, 


The dispersion is computed from (22), with account 


being taken of (24), (25), and (26), 





™m 
| ti+ > i*AyAg..-Ay_ 10 sie 
i=2 f y 
D{y| =a" ™m +2—,4 
1+ Sat DS) Aads...Ay_ ih 
1=2 


For m = 1, we easily get from (21) and (24) that 
M[y] = 1.5a?; the same result is also obtained for non- 











where, as before, Az and ©] are determined from (17) 


with account being taken of (I). 


We transform the first term in the braces. 
denote it by B(m), It is easily taken to the form 








2 3 i* ml2 
B(m) = = (m — i)\(m + i)! 
~ mi!2 
25) 





~ £4, (m — i)! (m + i)! 


The denominator is the same as in the analogous 
expression for M[y], andequals 2" ms? ; (2m)!, 


We now consider the series in the numerator. 


m ° 
\ . 
“J (m —i)'(m +i 


i=0 





i= mtA (m) — 


— (2m? — 2m -+- 1) A (m — 1) + A(m — 2). 


By substituting this expression in B(m) we get 


m (3m — 1) 


B (m) = 








It is obvious that all these formulas are valid for 


m > 1, 


Finally, for m > 1 ( n/4a* > 1), we find that 


D{y| = 


aS. 2 
63 2 ). 











For all n/ 4a” < 1, we immediately obtain from 


(22) and (24) that 


Dy] =- 


2. Mathematical Expectation and Disper- 











sion of the Output Variable for a Uniform 





Distribution Law in the Case of a Piece- 








wise-Linear Object Characteristic 





o[ 
For the case when the object characteristic is y =|x| , a 
the general expression for the mathematical expectation 


of the output is found in [1], We recall that 








aS if 
T=C+%, de (m0, 4.) tyeanl! 
3 i, 0 ( ) 1.0 10, if 
—ta 
(Av, |=2a; p= p= | pra 
—0o 
1 i-1 
n= i > fo = 2qr; rm, 
9 qrenmnenp 
e+ =) 
—1}—M? fy}, oti 
‘ 1 
Tv 447 — AS 
er M{yj=a a, 
q Pans on 
q rae 
In the case of the uniform distribution of (23) 
Then we a 1 a * 1 a (1 a \ 
jz-% (z>q) fat G>4) 
ies (Seth: Me 4 <4) 
(<5): \2 Sq)" 
a\2 a 1 
PY Fd oa a 
for — <a» 
M [y] émne2 4 2 
x rag 7 
a a 1 
=. fi = oe 
[ = or = > 7° 
We compute the dispersion 
{ 2 
—a a, a = — 
"Th... 
= et 
a ee 
y We 2 It 
D {yl = (Ay)r, — M*[y] = 2? (eee, 
oo el 
gr (i — A) +A 
=@ 1 2° 
— 4 qa 
(1— (q+ >) 
By substituting p and q, we get 
4s —0(Sy 
n i‘ a 1 
for —<=> 
a \2 n y ed 
D{y) = ti 
for £ >! .- 
" , 
EE 
Va for $ ee. 
0 for @~ 1 
“ ad 2 ‘i 


i=0, 
i +0; 
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THE ACTION OF RANDOM PROCESSES ON DISCONTINUOUS 


CONTROL SYSTEMS 


G. I. Pyatnitskii 
(Moscow) 


Translated from Avtomatika i Telemekhanika Vol, 21, No. 5, pp. 585-594, May, 1960 


Original article submitted April 13, 1959 


Questions of the statistical calculations of discontinuous (sampled-data) control systems are considered, 
these systems being acted upon by stationary random processes, We also consider one particular case when a 
nonstationary random process acts, It is shown that the method of statistical linearization can be used for the 
approximate calculation of nonlinear systems in the same form in which it is used for the calculation of contin- 


uous nonlinear systems, 


INTRODUCTION 


Today there is a heightened interest in sampled-data, 
or discontinuous, control systems, These include systems 
containing sampled-data elements of one form or another 
[1, 2], extrapolating devices [3, 4], and digital computers 
[5, 6). 

Discontinuous control systems are superior to con- 
tinuous ones in those cases when it is necessary to con- 
trol a slowly varying process, since they require less 
energy for their operation, 

Sampled-data systems can successfully be employed, 
and are employed, for the purpose of measurement, particu- 
larly telemetry; with this, they may easily be made 
multichannel devices, 

Sampled-data systems can be used for the control of 
several objects by means of the same apparatus which a 
continuous system must use to control just one object. 

The use of digital computers in control systems trans- 
forms them into discontinuous control systems since, by 
their nature, digital machines have to do with discrete 
information, with digits, Control systems with digital 
devices are very flexible since, by a simple change of 
program, they can be made capable of solving their prob- 
lem under different conditions, 

Discontinuous control systems of a quite broad class 
can be described from a single point of view. In this 
Class fall those systems whose behavior is determined by 
the values of the input signals at discrete equispaced mo- 
ments of time. 

In [7] a method is given for describing the operation 
of sampled-data elements with arbitrary pulse forms, The 
essence of the method is that the actual sampled-data ele- 
ment is presented in the form of a series connection of 
ideal sampled-data elements which process pulses in the 
form of Dirac 6 -functions (following one another at equal- 
ly spaced moments of time with norms equal to the values 
of the element’s input signal at the corresponding moments 
of time) and which form a link which transforms the 6 - 


pulses to pulses of the necessary forms, The shaping link 
is linear, and is described by a certain weight function and 
transfer function, 








oe ce ee ee 71 

| 

cnseeeneneupenmentinaneal 
Fig, 1. 


The series connection of the sampled-data element 
with a linear system forms a linear sampled-data circuit. 
The shaping link with the linear portion of the circuit 
forms the equivalent continuous portion (ECP) of the cir- 
cuit (Fig. 1), The transfer function of the ECP is, ob- 
viously, the product of the transfer functions of the shap- 
ing link and the linear continuous portion of the circuit, 
and the weight function of the ECP is the convolution of 
the weight functions of the shaping link and the continu~ 
ous portion, 

If x(t) is the signal acting on the sampled-data cir- 
cuit’s input y (t) is the output signal and g(t) is the weight 
function of the ECP then, in the steady state of the cir- 
cuit, the following relationship holds; : 


k 
y (kAt + eAt) = bi z(nAt) g ((k —n) At + eAt}. 
n=—oo 
(1) 
Here, the time t is given in the form t = kAt + €At, 
where k is an integer, € <1 and At is the interval of dis- 
creteness (quantization). 
Application of the discrete Laplace transform [1] to 
(1) gives 


Y"(p, ©) =G (p, 8) X°(p), (2) 
where Y* , G*, and X® are the transforms of y, g, x, re- 
spectively, For example, 

G*(p, ®) = >} g(mAt + edt) e—Pm', 
m=0 
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By analogy with continuous systems, G*(p, €) is 
called the transfer function of the sampled-data circuit. 
It may be shown that G*(p, €) is related to the transfer 
function (in the usual sense) of the ECP, G(p) by the re- 
lationship [8] 
on .2nm\ *4! p+i =e 
» G(p + ise ites 
If p is replaced by iw, one obtains the frequency 
characteristic of the sampled-data circuit, which may 
be interpreted as the ratio of the output signal to the in- 
put signal if the input signal is a function of the form 
elwnAt 
Relationships (1) and (2) also hold for systems with 
feedback, except that, instead of G* (p, € ) in (2), the 
system's transfer function enters and, instead of g(t) in 
(1), the properly defined weight function of system enters. 
In the sequel, a time-domain function will always 
be denoted by a lower-case Latin letter, while its trans~ 
form will be denoted by the corresponding upper-case 
letter; to distinguish discrete from continuous transforms, 
the discrete transforms will always be marked with a star, 


1. Stationary Random Processes with 
Discrete Arguments and Their Actions 
on Linear Sampled-Data Systems 











For the following exposition we shall need certain 
results from the theory of stationary (in the broad sense) 
random processes with discrete arguments, This theory 
is presented in great detail in the book of J. L. Doob [9]. 
In [10], the theory of discrete random processes is applied 
to a probabilistic computation of sampled-data systems, 
In contradistinction to [10], we shall here take as basic 
the spectral representation of stationary processes given 
in [9]. 

Each stationary random process x(nAt) with zero 
mathematical expectation, permits of a spectral repre - 
sentation 
n/At 

\ eion Aly (w), 


—n/At 


x (nAt) = ot 


z (4) 
where the random process u(w) has orthogonal increments, 
L.e., 


0 for 


| W; =F We; 
M [du (@,) du (2)| = {st aF (w) for @,; =@,=. 


(5) 
Here, M is the symbol for mathematical expecta - 
tion, and a superscript bar denotes the complex conjugate. 
The correlation function of a discrete process can be 
represented in the form 


n/Al 
R(ndt) = 34 \ elun Ald F (w), 


—nx/At 


(6) 





The function F(w) [the same function as in (5) ] is 
called the spectral function of the process, and its deriva - 
tive with respect to w, i.e., S* (w) = Fw), is called the 
spectral density. 

The spectral density can be presented in the form 
of a Fourier series whose coefficients are the values of 
the correlation function 

oo 
>; R (nAt) etonst, 


n=—0o 


5" (a) = (1) 
It follows from (7) that the spectral density of a dis- 


crete process is a periodic function of w with a period of 
am/ At 


A discrete random process may be obtained from 
a continuous one if one consider its values at discrete, 
equally spaced, moments of time. 

If x(t) is a continuous stationary random process, 
then its values for t = nAt comprise a stationary discrete 
process, Let S(w) be the spectral density of the continu- 
ous process, The spectral density S*(w) of the discrete 
process obtained from it can be computed by means of 
the relationship [11) 

. 1 
S (@) = Al 
™m 
which is easily obtained from (7). 

If the stationary random process x(nAt), which can be 
presented in the form of (4), acts on the input of some 
linear sampled-data system with constant parameters,then 
the system's steady-state output, i.e., after damping of 
the transient response in the system, can be computed by 
taking account of (1) and (4), The output process will 
have the spectral representation 


S(o a a ; 


> (8) 


co 


n/At 
\ ekAIG" (iw, ¢)du(w), 


—n/At 


At 


y (kAt + eAt) = = 
where G* (iw, €) is the system's frequency characteristic, 

It follows from (9), when (5) is taken into account, 
that the output process" correlation function is 


n/At 
At ; . 
Ry(ndt, edt) = 3 \ end’ (w) dw, (19) 
—n/At 
where S,*(w) is the spectral density of the process y, re- 


lated to the spectral density of process x(nAt) by the 
relationship 
S(o, ©) = |G" (iw, €) [*Sz(w). (11) 


The dispersion of the process y can obviously be com- 
puted from the formula 


n/At 
3 (e) = 3 \ |G (iw, €)[2S%(w)do. (2) 
—n/Al 


Formulas (10) and (12) were obtained from [10]. 
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We turn nowto (8), It immediately follows from it th 
that the ideal sampled-data element transforms the high - 
frequency region of the spectrum to a low-frequency one, 
This leads to the consequence that the introduction of a 
sampled-data element into a continuous system increases 
the random errors, This fact is seen particularly clearly 
in the following example, 

Let there be some continuous system, and let there 
act on it a stationary random process x(t) with spectral 
density S(w), Let the system's frequency characteristic 
K{iw) be such that for |w | > wo, its modulus is so small 
that the dispersion of the output process can be computed 
by the formula 


We 


a 1 . 
on = \ S(w)| K (iw) |? do. (13) 


— We 


Let us now connect this system to some sampled- 
data element whose pulse repetition period At satisfies 
the inequality 


At<—» 


The dispersion of the output random process from the 
sampled-data circuit thus obtained is computed from (12), 
With the assumptions made relative to K(iw) and At, the 
frequency characteristic of such a circuit will, according 
to (3), be 

° * = Ei 7 se . i eAt 

G* (iw, ¢) = Ar F (iw) K (iw) ef * 4) 

where F(iw) is the frequency characteristic of the sam - 
pled-data element's shaping link, 


We assume that over the entire frequency range where 
| K(iw) | still differs markedly from zero, | F(iw)| / At 


njAt 


differs little from unity. In fact, this assumption can be .» 
satisfied by choosing either a low-frequency part or a suf- 
ficiently small At, because | F(iw)| /At +1 for w0 
and differs greatly from unity only close to the frequencies 
w= +%/ At, 

Starting from this, we can write the following ex- 
pression for the dispersion at the sampled-data circuit's 
output: 

n/At co ° 

° mm 

|K(iw)*? S\ Sz (w " Fr) dw. 
—n/Ai m=—oo 
(15) 

By comparing (15) with (13), and taking into account 

that w, < */At, we can infer that of >on 


2. Action of a Nonstationary Signal on 
a Linear Sampled-Data System 








Let a nonstationary signal of the following form act 
on a linear sampled-data system 


E(t) =r (t)2(0), as) 
where x(t) is a stationary random process and r(t) is a 


function of time, 
The output process can be found by (1): 


k 
y(kAt+eAt)= >) &(ndt)g[(k—n) At + At}. 


The process € (nAt) can obviously be presented in 
the form 


n/At 
§(nAt) = r (nAt) \ endidy(w), (18) 
—n/At 


By substituting (18) in (17) we get 


y (kAt + eAt) = = \ du (w) ekAt S$) r ((k — m) At] g(mAt + dt) etomar, (18) 


—n/At m=0 


Expression (19) is the spectral representation of the output process. It is a nonstationary process, since k en- 
ters into the integrand, and not only into the exponential factor, 


We obtain the dispersion of the process y in the following form. 


n/At 
oy (k, &)= on 
—n/At 


We consider two particular cases: 
1, Let r(t) be a linear function of time: 


r(t)=a- bt. 


= \ daoS; («) | > r{(k — m) At] g(mAt + eAt) e-tomat . (20) 
m=0 






























Then, 
S) [a+ b(k—m) At} g [(m + 2) At] e-tomat — 


m=0 
= (a + DkAt)G* (iw, 2) + ib “ G (iw, 2). 
If the system's frequency characteristic G*°(lw, €) is 
presented in the form 


G’ (iw, ¢) = P (w) e*), 


then, for the square of the modulus of the sum under the 
integral sign in (20), we obtain the expression 


P? ((a + bkAt) 4- be’ |? + bP”, 
where the “prime® denotes differentiation with respect 
to Ww. 
We obtain the following final expression for the dis- 
persion: 


nr/At 
oy(k, &) = 5 dos; (w) (P? [(a + bkAL) + 
—n7/At 


+ bg’ + bP), 
(22) 


For computations by (22), it is necessary to differenti- 
ate the modulus and phase frequency characteristics with 
respect toW, These operations may be implemented 
graphically, 

2, Let 


r(t) = e**. (21a) 


In this case, 


summation sign outside of this neighborhood be so small 
that it would be possible to neglect it; at the same time, 
inside the neighborhood,r(t) must differ only slightly from 
a straight line or an exponential. 

By using general relationship (20), one can compute 
the dispersion in cases when r(t) varies by a more general 
law than (21) or (21a), 


3. Action of a Stationary Random Process 
on a Sampled-Data System which Contains 
a Nonlinear Element 











We assume that the stationary random process of the 
form 
A x/At 
t , » 
x(nAt + eAt) = = \ ebondt()* (w, ¢) du (w) 


—n/At (24) 


acts on some continuous linear system with weight func- 
tion k(t), The steady-state output process can be presented 


in the form , 


y (mAt +- yAt) — >) At\x [((m —n) At + 
m=0 0 
+ eAt) k[nAt -+- (4 —¢) At) ds. 
(25) 


By substituting (24) therein, we obtain the spectral 
representation of the output process 


m/At 
y (mAt + v,Al) = ~ \ du (m) e&™4t P* (@, 1), 
4-m/At 


co 

: 26 

m=0 1 
P*(m, 4) = A\ K*(@, 4-2) @'(, s)dz, @ 

and, consequently, 4 

n/At foe) 
oi, (k, ¢) — e@akAt = \ doS.. (w) |G" (a + iw, ©) 2, kK’ (@, a. e) = 2 k [mAt + (y, eared s) Ai| e—iomaAt 

~ ened — (28) 


(23) 

In order to compute the dispersion by this formula, 
it is necessary to reconstruct the system's frequency char- 
acteristic, since the presence of « would change all the 
time constants of the system's continuous part, 

A formula analogous to (23) for the case of continu- 
ous systems was obtained by V. S. Krapivin [12]. 

Computations by the formulas obtained can be car- 
ried out in the case when the function f(t) is neither linear 
nor exponential, For this, it is necessary that r(t), in the 
neighborhood of the moment of time ty of interest, can be 
approximately presented in the form of a linear or expo- 
nential function, The size of this neighborhood can be 
estimated on the basis of the sum entering into (20), It 
must be such that the magnitude of the quantity under the 
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We now consider the sampled-data circuit shown on 
Fig. 2, This circuit contains a nonlinear element, Such 
an element might be, for example, a relay element, an 
element with a limited linear zone, and other elements 
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Fig. 2. Nonlinear sampled-data circuit, 
1 is an ideal sampled-data element, 2 
is a linear link, 3 is a nonlinear element, 
4 is a linear link, 
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[13]. At the circuit's input acts a stationary random pro- 
cess v(t) = My + Vo(t), where my is the mathematical ex- 
pectation of the process v and vg (t) is a process with zero 
mathematical expectation. By using the results of section 
1, we can easily find the spectral representation and, con- 
sequently, the correlation function Rz (nAt, € At) of the 
process z at the nonlinear element's input. It is also easy 


to compute the mathematical expectation mz of the process 


z. We shall assume that the process vit) is normally 
distributed. In this case, the correlation function 

R, (n At, eAt) of the process x at the nonlinear element's 
output can be presented in the form of a series [14] 


co 
Rz(nAt, eAt) = 1“ pl (nt, edt). 
l=1 


(29) 


Here, p, = Rz / of, o% is the dispersion of process z, 
and a, are coefficients which depend on 97 and m,. In 
its turn, the mathematical expectation m, of process x 
can be presented in the form 7 


m,y=y¥ (S25 m;) Mz. 


(30) 
A simple and elegant method of calculating the co- 


efficients al and x, based on the use of the 4 -function, is 
presented in [14], 


For an element with a limited linear zone, 
b { Sc, a—_ m, a—m 
E+ BE 0 
% a m, z cS, 


a-+-m 


5) 9a) +2ra) 


os) — ot (SP. 


“5 




















OF Ft 


where a is the magnitude of the element's linearity zone, 
b is the limitation level and b/ a is the slope of the linear 
part of the element's characteristic, 


ps — 


For a relay element, 


x= 5,5 [20 (=) —4], ! 
a, = 4b? [ o” (— “yf 


where b is the signal level at the relay element's output. 
x ig 


1 _ 
® (x) = —— é dt, 
V2n 2. 
and @(1) (x) is the [th derivative of (x) with respect to 
x; 9 (2) = ®' (2). : 
The spectral density of process x, starting from (29), 
can obviously be given in the form of a series 


Here, everywhere 


Sx (w, ©) = a,S;(, &) + 34 ar zara So (w, ®), 


(31) 
where So (u, €) is the spectral density of process z, and 
S? (2) is the ? th-order spectrum of the function se (w, €). 

Tables 1 and 2 give the values of the coefficients 
a/l! for an element with a limited linearity zone and 
a relay element for a = b = 1, for various values of m, 
and 9,, 


It is clear from the tables that in the overwhelming 
majority of the cases there holds the inequality 
a, /U! << ay,for! = 2,3,... Thus, the terms with iz 2 
in (31) can be ignored, 

We turn back now to Fig, 2, Let link 2 have weight 
function g(t) and link 4, weight function k(t), Then the 
spectral density of process z will be 

S3(w, &) = |G" (iw, e)|? 5; (@). (32) 

As a consequence of what has already been said about 

sy (iw, €), we have that 


Sz (o, &) + a,|G" (iw, ©) 25; (@). (33) 








TABLE 1 

m, Sc, a, a,/2 a,/6 aJf%& 
0 0.5 0.2220 0 0.00194 0 

0 1.0 0.4640 0 0.03900 0 

0 2.0 0.5860 0 0.08270 0 

0.5 0.5 0.1760 0.0070 0.00272 0.00001 
0.5 1.0 0.3900 0.0248 0.02280 0.00754 
0.5 2.0 0.5540 0.0118 0.06950 0.00916 
1.0 0.5 0.0625 0.0198 0 0.00167 
1.0 1.0 0.2270 0.0594 0.00194 0.01310 
1.0 2.0 0.4650 0.0494 0.03900 0.02650 
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TABLE 2 




















m, 6, a, a,/2 a,/6 a,/24 
0 0.5 0.6360 0 0.1060 0 

0 1.0 0.6360 0 0. 1060 0 

0 2.0 0.6360 0 0.1060 0 

0.5 0.5 0.2340 0.1170 0 0.03910 
0.5 1.0 0.4950 0.0620 0.0495 0.03910 
0.5 2.0 0.5970 0.0177 0.0865 0.01340 
1.0 0.5 0.0116 0.0233 0.0175 0.00194 
1.0 1.0 0.2340 0.1170 0 0.03910 
1.0 2.0 0.4950 0.0620 0.0465 0.03910 














From this we find that the approximate spectral repre- 
sentation of the process x has the form 


n/At 


x (nt +- cAt) = = \ cont VY a. G" (iw, ©) du (w), 
—n/At 
(34) 
. 1 as — ae os . 27 
P* (@, 4) = = » V a,K (iw 7t = 


This expression is the frequency characteristic of the 
circuit in the case when the nonlinear element is replaced 
by an amplifier with gain equal to Ya, 

By using the representation of the nonlinear element 
in the form of a linear amplifier, we can easily calculate 
the probabilistic characteristics of the system with feed- 
back, It is assumed for this that the process at the non- 
linear element's input is normal. Such: a method of cal- 
culating continuous systems is known as the method of 
statistical linearization, developed by L E. Kazakov[13]}. 

Thus, the statistical linearization method is also ap- 
plicable to sampled-data systems, The computational 
methods for the discrete case are no different from those 
used in the continuous case, which are well-known. The 
coefficient 4a ,is computed for different types of non- 
linear elements in [13], 


CONCLUSIONS 


In this paper we considered a number of questions 
dealing with the statistical computations of discontinuous 
control systems, We took as fundamental the spectral : 
representation of stationary random processes with dis- 
crete arguments, 

Starting from the relationship between the spectrum 
of a continuous process and that of the discrete process 
related to it, we found that under the same conditions the 
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)G (iw + i) aie ta 


where u(w) is the spectral process of vg(nAt). 

Process x acts on the input of linear link 4, We find 
the spectral representation of process y at the circuits out- 
put in the form of (26). 

By taking into account that, in the given case, 


Q° (iw, ¢) 8 Ya, G (iw, ©), we obtain 


2nn ) (35) 
Al 


dispersion at the output will be larger in the sampled-data 
system than in the continuous system equivalent to it, 
This is a consequence of the fact that the sampled-data 
element transforms the high-frequency part of the input 
process spectrum into a low-frequency region, 

The formulas obtained for the dispersion at the out- 
put of stationary systems permit the computations to be 
carried out in the cases when nonstationary processes act, 

An investigation of the question of the action of a 
stationary process on a discrete system containing a non- 
linear element led to the conclusion that the statistical 
linearization method is applicable to the discrete case, 
the computational methodology remaining the same as 
in the case of continuous systems, 
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ERROR ACCUMULATION IN DIGITAL COMPUTATIONS 


V. A, Brik 
(Moscow) 


Translated from Avtomatika { Telemekhanika Vol, 21, No. 5, pp. 595-600, May, 1960 


Original article submitted December 26, 1959 


The paper describes a method, based on diffusion theory, of estimating the errors of computing devices 
used in automatic control systems, as well as special-purpose and general-purpose digital computers, An investi - 
gation is made of the errors of cyclical computational processes, in which the successively accumulated ele - 
mentary errors can be considered as independent random variables with an approximately normal distribution 


density. 


Digital computing devices which solve recursive dif- 
ference equations [1] are used in automatic control sys- 
tems, In the operation of these devices, the computations 
are implemented cyclically, wherein the results of the 
previous cycle, or cycles, are the input data for the follow- 
ing cycle, For this, it is desirable to have the results ob- 
tained after each cycle, Recursive computational pro- 
cesses with output of intermediate data also arise frequent- 
ly in the solution of various problems on general-purpose 
or special- purpose computers, for example, in computing 
tables of some function or other by recursion formulas, 

With each computational cycle there ordinarily arise 
errors which gradually accumulate, This accumulation 
can occur either quite regularly or randomly, If the suc- 
cessively arising errors are virtually statistically independ- 
ent, then they may be estimated by the use of probabilistic 
methods [2, 3]. The possibility of using these methods can 
be established ['3) after an analysis of the statistical prop- 
erties of the table of successive errors. 

In many cases, the process of error accumulation 
amounts to the addition, to some total error which de- 
termines the quality of the entire solution and which exists 
prior to the given cycle, of an error resulting from the 
given cycle, 

If the error attributable to the given cycle is the re- 
sult of the effect of several independent sources of error 
(for example, rounding errors after several multiplications), 
then the resulting error of the cycle may have a probabil- 
ity density close to normal, Ordinarily, the mean value 
of the error in .a cycle is zero, The dispersion 9%? of the 
cycle's error is computed from an analysis of the table of 
successive errors of the cycles, 

If one is given some limit a > 0 which must not be 
exceeded by the absolute magnitude of the total error 
4. then, if the conditions stated above are met (addition 
of. the successive errors of the cycles, statistical independ- 
ence of these errors), the process of error accumulation 
can be considered as the random walk (diffusion) of a 
point between two absorbing barriers (Fig. 1). The ordinal 
numbers of the successive cycles are laid out along the 
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Fig. 1. Diffusion process 
with two absorbing barriers, 


axis of abscissas on Fig. 1. The total error after the nth 
cycle is the sum of the individual cycle errors, 


An = >) &. 
i=1 


The random walk begins at the point A, = 0 (it is 
assumed that the initial data are given exactly). We shall 
assume that all the increments 4, are independent ran- 
dorn variables with identical normal distributions with dis- 
persion 9” and zero mean value. 

We now find the probability P,(n) that, in n cycles, 
the random walk will go beyond the barriers at least once’ 

For a given n, the probability density of the random 
variable Ap has the form 


f(a, n) = rove? (— =). @ 


Using the reflection method [4, 6], we can introduce 
a modified "distribution density” 


fle m) = (A aig one ee 


s=—0O 


(2) 
characterizing the probability distribution law of the com- 
plex event 


f(z, n)dx= P{xQAngz+ dz; 


| Ai] <a} (i =1,2,.... n—4). 


* For the case of continuous n, the problem,in a more 
general form,was solved by Bartlett [4, 5]. 
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With this, consequently, 


| f(z, n)dz = P {| An| <a: 


|Ai| << a} (i =4, 2,..., 2 —1). 


| fle, n)dz = P{|Ai| <a) 


—a 


(é = 1, 2, ..., a). 


The probability of interest to us, Pg(n), equals 
Py(n) =1— fi(x, n)dz. (3) 
—~a 


After some transformations, we can obtain from (2) 
and (3) that 


Pyin)=45) Fi} (1 = 


=4{F(=2)— F (=) 4 F(=*)-..}, (4) 
where 


u 


F(u) = ( 7me (— +) dz. 


—co 

Po(n), as a function of a/© and n, computed by (4), 
is shown on Fig, 2, 

The method of estimating the error consists in find- 
ing such a value of a that, for a given number of cycles 
n and a dispersion o® of the maximum possible error 
Dn (Dp is the maximum of all the Aj, i = 1,2,...,m), 
Pon) will equal some definite, sufficiently small, quant- 
ity a (for example, a = 107%, 107° or 10), Thus, a 
is the confidence level for which 


—Da< hic Dn (i =4, 2, ..., n). 
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The dependence of D,/® on n for different a is 
shown in Fig. 3 on a semilogarithmic scale, 


The graphs of Fig. 3 answer the question as to the 
practical limits of the errors A; (i= 1,2,...,Mm) fora 
known dispersion 9? of the elementary errors (cycle er- 
rors), These graphs can be used for estimating the errors 
of any computing process in which there occurs accumula - 
tion (addition) of statistically independent successive er- 
rors with mean zero, dispersion 97, and a distribution close 
to a normal, 


We now make an estimate of the error of computing 
the coordinates of successive equally spaced points on a 


Dn/o 


ry 





Fig. 3, Dependence of the maxi- 
mum error D,, divided by the 
mean-square deviation of the in- 
crement 9, on the number of cy- 
cles n and on the quantity a, 





=i 


120 0 160 180 200 a/8 


Fig. 2, The probability Pgn) of going beyond the barriers as a func- 
tion of the number of cycles n and of the ratio of the distance to the 
barriers a to the mean-square deviation of the increment ©, 
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circle (with center at the origin) by means of the recursion 
formulas 


x(n +1) =a(n)-(1— 2°) Fy (n).2™, 
y(n + 1) =y(n)-(1—2-*) + a(n).2, 


which are obtained from the formulas 
z(n + 1) = 2z(n)cosw+ y(n) sina, 
y(n + 1) = y(n) cosw+2(n) sino 
after the substitutions 


coso~1— 27"! sinw = 2", 


(5) 


where k is a positive integer and w is a constant angular 
increment, This algorithm is used in a special-purpose 
computer [7] which prepares programs for automatic mill- 
ing machines, 


The error to be investigated will be the deviation of 
the current radius Ry = V 7° (i) + y*(i) from the initial 
radius J%, = VY x?(V) + y®? (G6) over the entire circum - 
ference of the circle. We shall only investigate the errors 
AR; = Rj — Rg arising from rounding after multiplication 
(in binary form) by the quantities 272-1 and 27k, The 
methodical error arising from the substitution in (5), for 
sufficiently large k (larger than three to five), is small and 
may be neglected, 

We shall denote the rounding errors arising from the 
multiplications x(n)27k~1, y(ny2-k, y(ny2-**"1, and 
x(n)2-k by € ,(n), € o(n), € x(n), and € ,(n),respectively, 
where an error leading to an increase of the radius will 
be considered to be positive. 

It can be shown that the error of a cycle equals 


4 = Ri — Re, = li, (i— 1) + & (§— 1)] cos 8 + 


~ [is (i — 1) + Se - 








1)} sin 3, 





- —S TFS 





where 


i i lee 
3 = AFC, 

We shall assume that the computer operates in fixed 
point, that !z|< 1, |y!< 1, and that the machine uses 
an N-bit word after the binary point. A statistical analysis 
of the successive errors € , carried out for the given prob- 
lem, showed that for A> 2-\*4**?  ive,, for sufficient- 
ly large N, all the quantities € ,(n), & ,(n), & s(n), &4(m) are 
statistically independent, and successive values do not de- 
pend on preceding ones, It turned out, with this, that they 
have a uniform distribution density P(2) (Fis. 4). The dis- 
persion of such a distribution equals 2” /12, and the dis- 
persion of the random variables 4 j = Ri ~Rj-y. the sum of 
the four independent random variables, equals 


9—2N 9—2N —2N 


= — 2 ) <- . 2 _ 2 
o% =22 77 098" B + 2 5 sin p=--. 





The graphs 


of the distribution density p(z) of the quantity 5; for vari- 
ous 8, computed by a convolution integral, are shown on 
Fig. 5. Also constructed on the figure is the density curve 
for a normal distribution with the same dispersion, This 
curve differs slightly from the others, so that one can con- 
sider that all the 5; are normally distributed with disper- 
sion 27° /g, 


p(2) 


2 
z 

















-2°4-1 0 26! Z 
Fig. 4. Distribution density of 


the quantity & (n), 





~2 4 0 ‘ 
fig. 5. Distribution density of the increment 5; for various 8: 1 is for 8 = 
ma /2; 2 is forB =m™/2 + 7/4; 3 is for8 =mn/2+ 1/8, m= 0,1, 2, 3; 


4 is a normal distribution, 
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The number of cycles n for the compiete circle, 
k+1 f f 
equal approximately to 72""", changes by a factor o 
two when k changes by unity (cf., the table), 
By using the graphs given on Fig, 3, one can deter- 
mine the dependence of D, onk, For example, for k = 7 
and for a selected a =10~* the maximum error 


i 

| 42% | max= | >) 8] = D, (i = 1, a cee n) 

j=) | Max 

equals approximately (40)2-N, tt is of interest to com- 
pare this value with the (theoretically) maximum possi- 
ble error AR,. Since the maximum value of Ri -Rj-;, 
for different 8, varies from 27N to hb S aed and, 
on the average, equals (1,3)2 “No then the theoretically 
possible error AR, equals 1.3n27 = (1040)2-N, which 
is 26 times the value obtained for Dp. 

Instead of using graphs to determine Dy, one can use 
the well-known tables of the function y = F(u). By sub- 
stituting « instead of Pan) and Dy insteadof a in(4), 
and by discarding all except the first term, we get the 
equation 





a= AF ( me) . 
o Vn, 
Knowing &, o, and n, we can easily determine D, 
from the tables, 
The simplification of (4) just given is admissible for 
D,/¢ Vn, ‘greater than 1 to 1.5 since, with this, the 
sum of the discarded terms, as a fraction of the first term, 


does not exceed 0,01 to 0.00005, For small Pg(n) (or «), 
this condition is always met, 

In about the same way, one can estimate the accumu-~ 
lated errors in those cases when the normal distribution of 
the independent variables 5, has a nonzero mean. 
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- 


The effect of velocity feedback on the free oscillations of a loaded relay servomechanism is considered. 
The servomechanism's linear portion is described by a complete second-order equation. The investigation is 
carried out by the point method, using the theory of point transformations [1]. 

The dependence of the amplitude and time of the autooscillations* half periods on the system's parameters, 


including the velocity feedback coefficient, is determined, 


1, Introduction and Posing of the Problem 





The investigation of the dynamics of relay servo- 
mechanisms whose linear parts are described by complete 
second-order equations, has been carried out until now 
for the cases when the servomechanism’s controlling ele - 
ment has a symmetric nonlinear characteristic, In [2-5] 
there were considered servomechanisms with Z-shaped or 


looped characteristics and, in [6], characteristics with loops 


and insensitive zones, 


Few works have been published [7] in which considera - 


tion is given to unsymmertric looped characteristics of 
loaded-servomechanism controlling elements, 

Some automatic control systems use electrical servo- 
mechanisms in which, by means of special relay circuits, 
there is implemented a charge in direction of the veloc - 
ity by the switching of the supply voltage of an independ- 
ently-excited motor from 24 to 48 v. To decrease the 


static error due to static friction, the servomechanisms 
are implemented in the form of autooscillatory systems, 
The autooscillations arise due to the unsymmetric char- 
acteristic of the controlling relay, which has one normal- 
ly closed contact, In the servomechanism whose schematic 
is shown in Fig, 1, both direct and velocity (position and 
rate) feedback are employed. 


2. Equations of Motion 





As is clear from Fig. 1, the total voltage applied to 
the circuit is 


AU = Up, — Upp — Use, (1) 
where Ujn is the input signal voltage, Up, is the (direct) 
feedback voltage and Uys is the velocity feedback voltage 
(from the tachometer generator 4), 

Controlling relay Ry, connected in the tube's anode 
circuit, has two working contacts 1R,; and 2R;. With no 
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current in relay winding 1R, (with negative grid bias), 
the relay closes contact 1R by means of a spring. With 

a positive grid bias, anode current will flow through the 
tube, causing relay R, to operate and closing contact 2R;. 
In this case, when contact 1R, closes, current relays Ry 
and R, and distributor relay R, are simultaneously in ac- 
tion, If contact 2R, closes, then current relays Rg and Ry 
and distributor relay R, will be simultaneously in action, 
It is necessary to mention that, with zero signal on the 
input, relay R, operates with a frequency of 20 cps, With 
this, relay R, closes its working contact 1R, with the dou- 
bled frequency 40 cps, since it is switched-in twice as 
frequently as the current relays, When relay R, operates, 
its armature touches contact 1R, and switches in coil 1R, 
of current relay Rg. Operation of current relays R, and R, 
leads to a rotation of motor 1 in one direction, while 
operation of relays Rg and R, cause it to turn in the other 
direction, By its operation, relay Rg switches the motor 
armature from the 24 v source to the 48 v source, When 
the winding of relay 1R, is disconnected, its contact 1Rg 
closes the motor armature circuit on that side of resistor 
rt which is connected to part of the 24 v storage battery. 


When relay Rg operates, its contact 2R, closes the 
motor armature circuit to the other part of resistor r, 
which is connected to the full voltage of the 48 v storage 
battery, With no input signal, the motor will still be ro- 
tated by virtue of the absence of a neutral position for re- 
lay Ry. With zero input signal, the process of switching 
relay R, occurs so rapidly that the contact portionof re- 
lay Rg, due to its own inertia, cannot operate, and relay 
Rg remains unconnected, With this, current relays Ry, Rs 
or Rg,Ry Operate, feeding 24 v to the motor armature, In 
this case, steering organ 2, set into motion by the servo- 
mechanism, will oscillate with a frequency of 20 cps and 
an amplitude of about 0,5°, When a command is applied, 
the contact of relay 2R, is closed. In this case, the dura- 
tion of the closure is sufficient to allow relay R,, and then 
relay Rg, to operate, In operating, relay Rs applies 48 v 
to the motor armature, With this, the motor will execute 
the command with maximal speed and torque, Thanks to 
the presence in the circuit of the tachometer generator 
feedback, the operation of relay R; will be accelerated, 
With this, the motor will be switched to 24 volts upon 
reaching the position given by the input signal. This will 
preclude the arising of a lengthy oscillatory process, 


We now consider the servomechanism's free oscilla - 
tions, i.e., its motion for Ujn = 0. The error equation of 
(1) can be given in the form 


AU = — Kt — K2Q,, (2) 


where « is the angle of rotation of the reducer*s output 
axis, 2, is the motor's angular rotational velocity, Ky 
and Ko are proportionality factors, 


The speed of the reducer'*s output shaft is related to 
the motor velocity by the relationship 2 = Q,/ K, =da /dt, 
where K, is the reduction coefficient, 


The equation of the electronic amplifier can be given 
in the form 


Ua= KaAU, (3) 


where U, is the voltage at the amplifier's output and K, 
is the gain, 


The equation of the relay block has the form 


U =P (Ug) = |U| F (Ua), (4) 


where U is the voltage applied to the motor circuit. 

We take the motor's voltage supply to be constant 
and equal to 48 v, i.e., we do not take into account the 
action of the scheme for switching the supply voltage. 
Such an assumption significantly simplifies the solution 
of the problem without leading to any essential errors in 
the investigation, since it is based on the fact that switch- 
ing of the supply voltage to 24 v occurs close to zero 
error, and turns out to he essential only for the autooscilla- 
tory mode, Consequently, in calculating the autooscilla- 
tion parameters, it is necessary to consider the motor’s 
supply voltage as equal to 24 v, while for the investigation 
of the system's transient motions one should take the sup- 
ply voltage as equally 48 v. 

The nonlinear function F(U,) is shown in Fig, 2, where 


(e + >) and (e — >) are the operate and release 
voltages,respectively, for the relays. 

F(l,) 

oubef 


a 








| 1 0 
TeH| Peay 


= ~--f-1 
Fig. 2, 


The equation for the motor's electrical circuit, if 
one neglects the armature’s inductance, looks as follows: 














U=IR+C,Q,, (5) 
where I is the armature current, R is the active impedance 
of the armature winding, E 7= C, 2, is the motor's counter- 
emf and C, is the motor constant, 

If we take the load torque M, to be proportional to 
the angle of rotation of the reducer’s output axis a, i.e., 
My = K,&, then the equation of motor motion can be 
given in the form 
dQ 
Jre a= 
where Cy is the motor constant and J,¢ is the moment of 
inertia, referred to the reducer’s output axis. By elimina- 
ting from (2)-(6) all the variables except U, = K,Ki% 
we get 


CyIK, — Kye, (6) 


d*Uu aU . . dU 
AT + BS} +CU,=—F(U,+D52),M 
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where 








, JR - C.Kr 
~ UTC, Kr Kak, * ~~ ~ [UTK,K,’ 
RK, KK 
Cm L , Dan o’r p 
|\U|C A,X, Ka Ky 


We introduce the following dimensionless param - 
eters; 


° C B 
t= ’ = ’ B— —_— , 
Vighe= Uy h 2V AC 


k=DY 5 y=C,d=0. 8 @) 


Equation (7), written in dimensionless quantities, has 
the form 


fat 2h 4 2——O(2+kH). 


The function @(°), where © = x +kdx /dt®*, is shown 
on Fig, 3, This function is described by the equation 





1 for ¢>u—-+> if @(o,)=1, 
®(c) = A 
— Lot o<4t sz if @(o) = —41, 
(10) 

where ©, is the initial condition at time t = +0, 

a) 

dra = 

= — A —e 




















Figure 4 shows the block schematic of the servo- 
mechanism scheme described by (9). Thus, the system 
is characterized by four essential parameters h, k, H¢, 


and A, 
-C6 "| fol 1 | 
p*+2hp+! 


Ir 
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Kp 
Fig, 4, 


3, The Phase Plane 





By setting dx/dt* = x= y, we write the equation of 
motion (9) in the form 


z=y, yt2hy+c=—O(r+hky). (11) 


We now consider the structure of the servomechanism's 
phase space, It follows, from the system of equations in 
(11), that the phase space of the system being investigated 
is a two-sheeted phase plane, On sheet I the system's 
motion is described by the differential equations 


zt=y, yt+2hy+2=—1. 
On sheet @, 


(12) 


t=y, y+2hy+r2=1. (13) 
From the system of equations (12) and (13), it follows 
that the system's equilibrium states are found at the points 
x=-1,y=0and x= 1, y= 0, which are stable foci (if 
h < 1) on stable nodes (if h > 1), Since, in actual systems, 





Ye 5 | <1 and u> 2 , then the representa- 


tive point, moving along a sheet of the phase plane, al- 
ways reaches the switching line, never falling on a point 
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of rest, For the purpose of system stabilization, the damp- where 

ing coefficient k is ordinarily chosen to be positive, The 

motion of the representative point along the two-shetted ' ( A ( A 
t,)=w —)— > 

phase plane occurs in the following manner. While mov- /(h)= wilt +H, : +h +> ) 

ing along sheet J [D(x + ky) = 1] (Fig, 5), the 

representative point, during some interval of time, falls 


on the line L(o=me+ +) from the line 


Ly (2 = — 7) a(t) = —w,(1+u,+4)+(1+n,—4) 


(w, cos w,t, aa Asin w,t;) emhts ’ 


The system's phase space is not symmetric, since its 
equilibrium states are found at points on the axis y = 0 
which are unsymmetrically placed with respect to the 
pointx= He, y= 90, 


(w, cos w,t, on h sin wyt,) eft, ’ 


d®y _ (k — 2h + 4/k) fs (ts) sin*wyty 
du? #3 (ts) ' 





4, Point Transformation (15) 





The problem of seeking limiting cycles and of in- where 
vestigating the partitioning of the phase space leads to the _ A*)., 
study of the point transformation of line L into itself, I(t) = — <> or [C(t +t 4 = |sin hs “Tr 
This point transformation T is defined ag T = S*S~. 
Transformation §* corresponds to a transition of the A\? A\? 
representative point, lying on line L, to a point lying on 7 (1 as ali 5) arr (! T Be + +) Pr 
line Ly, and transformation S~ corresponds to a transition 
of the representative point lying on line L, to a point 
lying on line L, We shall find the analytic expressions 
for the transformations St and S~, 


If the initial conditions are changed on sheet I, i.e., if 
X and y» are changed, the representative point will be 
translated along sheet I during different time irtervals, 


y 
Sheet I of the Phase Plane / [@M(c) = — {4} 








The case of small damping (h < 1). If, at time 
t® = 0, the representative point begins its motion from 
point (Xo, yg) of line L (cf, Fig. 5) and, during the time 
interval t® = t,, arrives at point (x,;, y,) of line L,, then 
the analytic expression for transformation S*, in the form 
of a correspondence function u, = u,(v), in the case when 
h < 1, will have the form 


A A 
(1 +p.+ 3 Jews - +¥.— | eltts 





























— w,-++ 
(k—2n+ >) \ sin wily 
>) 
(d+, + 7) (h—4&) A 
a a > 
k—2h+-> Fig. 6. 
14 
“ The transition time t, for the movement of the repre- 
—(144,-F) cos with + (1+ 4_+ +5 ea hts sentative point from line L to line L, by a trajectory ly- 
r -. ; w+ ing on sheet I and not going off on line L, i.e., the transi- 
( k—2h + *) sin wt, tion time along the "cylindrical® surface, lies in the inter- 
val (0, tys), The time T,, is defined by the condition that 
: t d 
(14,-4)@—0) trajectory m and line L be tangent at the point M(x,,, 
£2 +a— A Yav ) (cf., Fig. 6) 
1 ’ ya 
k— 2h a k : . 


ws (1+u—s)—(1+4, +4) 
where wi = 1-h?, u, = x», v = xq, and t, is a parameter, wert ? 
It follows from (14) that —At; 
(w, cos wit, + hsinw,ts)e *“=0. (16) 
do _ fi(tr) 


du, — fa (ts) ” The coordinates xagy and y,y of point M are also 
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defined by the condition that trajectory m and line L be 
tangent: 





(ssn) 
ites) 
lay k i - (A. + 
k—2h+ - 
(17) 
A 
1 r BP, + y 
OF ’ Yay ang x 4 . 
k—2h + = 


The transition time along a trajectory in sheet I and 
going off the L varies continuously from ty; to ™/w, as 
Xo varies, An investigation of the function v = v(u,) for 
the case when h <1 is given in Appendix I. The curve v = 
= v(uy), given by the parametric equations (14) and (15), 
is defined for all t, from 0 to ™ /wy. It is clear, from an 
analysis of the phase plane, that from the entire set of 
arms of the curve there should be retained only those arms 
which satisfy the following limits of coordinate variation: 


The curves v = v(u;) for values of system parameters 


such that h< 1|}y.+ + | <1 are given in Fig, 7. 





The case of large damping (h> 1), The transforma- 
tion S* in the form of a correspondence function u, = u,(v), 
for the case when h > 1, has an analytic expression which 
may be obtained from (14) by replacing wy, sin w,t, and 
COS Wyt; by, respectively, w, sinh wt, and cosh wt;, where 
w? =h?-1, In the case when h > 1, the transition time 
along a trajectory in sheet I which goes off line L varies 
within the interval (tys, a0). The time tys, in the case 
when h > 1, has an expression which may be obtained 
from (16) by replacing wy, sin wyty and cos wy,t, by, respec- 
tively, w sinh wt; and cosh wt;. For the case when h > 1, 
the coordinates x,, and y,,, are expressed, as in the 
case forh < 1, by (17). An investigation of the function 
v = v(u,) for the case when h > 1 is given in Appendix I, 
The limits of variation of u, and v are the same as in the 
case when h < 1, The curves of v = v(u,) for values of the 
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Sheet II of the Phase Plane [(@(c) = 1] 





The case of small damping (h < 1), If the representa - 
tive point, at time t®* = 0, begins its motion from point 
(xy, Y) on line L, (Fig. 5) and, after a time interval t®* =t,, 
arrives at the point (xg, y2) of line L, then the transforma - 
tion S~ in the form of a correspondence function ug, =u,(v) 
in the case h < 1 has an analytic expression which can 
be obtained from (14) by replacing v, uy, ty and H ¢ by, 
respectively, —Ug, —v, tz and ~ He, where uy = Xg, V = xy 
and t, is a parameter, 

The transition time t, for the movement of the repre - 
sentative point from line L, to line L along a trajectory 
in sheet II which does not leave line L, lies in the inter- 
val (0, tgs), The transition time by a trajectory in sheet 
II which leaves line L; varies continuously from tgs to 
 / wy as x, varies, The time tgs, defined by the condi- 
tion that trajectory n and line L, be tangent at the point 
N(Xgay, Y2 ay)» has an expression, in the case when h <1, 
which can be obtained from (16) by replacing H ¢ and ty, 
by, respectively, -H~ and tg, The coordinates x,,, and 
Y2 ay im the case h < 1 and in the caseh > lhave ex- 
pressions which can be obtained from (17) by replacing 
Xav » Yay » and HB ¢ by, respectively, —xgay. —yeayand 
“He 

“Investigation of the correspondence function u, = 
= u,(v) in the case h < 1 can be carried out in analogy 
with that of the function v = v(u,) in the case when h < 1, 
which is given in Appendix, I, by replacing v, uy ,Pe, ty 
tys and ty, by, respectively, —u,, —v, “HC + tap ty, and typ, 





From an analysis of the phase plane (Fig. 5), it is 
clear that of the entire set of arms of the curve u,=u,(v), 
there should remain only those arms which satisfy the fol- 
lowing limits of coordinate variation: 


Ug from Xgay to — 0, V from Xyay to w. 


Curves of u, =u,(v) for system parameter values 
h<i, 
shown on Fig, 9 are curves u, = u,(v) for system parameter 
values h < 1, |. +3 <1. 


The case of large damping (h > 1), The transforma- 
tion S~ in the form of a correspondence function ug= ug(v) 
in the case h > 1 has an analytic expression which can be 
obtained from (14) by replacing v, uy, €, Wy, ty, sin 
Wy,t, and cos wyt, by, respectively, —u,, -v, ~H €, W.te, 
sinh wt, and cosh wt,, In the case h > 1, the time t,, 
has an expression which can be obtained from (16) by re- 
placing He, wy, ty., sim syty, and Cos Wytys by, respective- 
ly, —Be, W, tgs, sinh wtgs and cosh wtp. 

Investigation of the function u, = u,(v) in the case 
h > 1 can be carried out in analogy to the investigation, 
given in Appendix II, of the function v = v(u,) in the case 
h > 1 by replacing v, uy, H ¢, ty, tys, tyg and ty by, re- 
spectively, —tlg, ~V, ~He, tg, tgs, tes and ty. The limits 
of variation of u, and v are the same as in the case h < 1, 
The curves ug = u,(v) are shown on Fig, 10 for system 


wet |<1. 





ee+ 3 | <1 are shown on Fig, 9. Also 








parameter values h > 1, 


h close to zero 
u,,Uy u,,U 


Uy, 











h close to uni 
Uy,U Y 


U,,Uy hb 2 

















Uy,Uy 
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Figure 10 also shows curves of u,; = u,(v) for system parameter values h > 1, 





A 
en 2 =a < . 
System Autooscillations 





It follows, from the diagrams of point transformation shown in Figs, 9 and 10 for the parameter values h == 1, 


pet =| <1, that for these system parameters, the curves of the correspondence function u, = u,(v) and u, = 





= u,(v) must always have one point of intersection in the first quadrant. Consequently, there exists one stable limit- 
ing cycle, i.e., there are always autooscillations in the system, 

These autooscillations can be established with the presence of a skipping mode (cf., Fig. 6). On the "cylin- 
drical surface” between lines L and L;, points MNyNM, mark out the region in which a skipping mode occurs, From 
the conditions that the curves u, = u,(v) and uy = u,(v) intersect, in the case h < 1, we find the times tig and tea of 
the half periods of autooscillation in the case when h < 1, These conditions are written in the form 


A\ ~ht 4a 
(t+u.+-)e (144, —+) cos witia 
sin w 


ee is oi i—» + cOS W, fy, 
2 2 


. 1» 
sin w, to, 





w, + 2(h—k)= 
14a 





(18) 
A A\ ht 
(1 +P, +) 008, ty, — (1 +u—sZ)e la 
sin Wi bia 





w, + 2(h—k)= 


A A)\ —nt 
(1—», —) c08 1, t59 — (1—», 4 Se = 


Ww, 
sin w, toa 





The amplitude of velocity autooscillation y is expressed, in the case when h < 1, by the equation y, = tt or 











4 (1+ me +) cosw,t,,—(1 + 1, —) eis 
“spage os “— 
= A - _A htgg 


—(1-1,+4)@—m], ™ 


where ty, and tg, are defined by (18), Analogously, for the case h > 1, the times of the autooscillation half periods 
tia and tg, and the autooscillation amplitude are found from the equations 


(1 +p. + a e ie _ (1 +t =) chwt,, 




















shwt,. w+2(kh—k)= 
A\ ht A 20 
(1—»,—F)e (14,4 >)chuty, = 
= w, 
sh wi,, 
A-wik<hew ‘ 
Uy ‘ P 
O<k<h-w k>hew } Wa 
| 
! 
| 
| 
! | 
1 
' 
22 aw v 
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A A 
(1 +eRe+ +) chwt,, _— (4 +,—3) ee 





sh Whig 


ti (1—»,-}) chur, —(1—1, + Fe 


—hloq 





sh wte, 





uy 


A A\ » 
5) ch wt, — (1 +u—p)e - 





{ [= 





¥, = 2k (k— 2h 44 sh wt,, = 
A a oes 
“i+n+$)n—y—lonrdese on ta 
—(1-1,+ 4) a=W |. 
SUMMARY 


With a looped relay characteristic of its controlling element, a loaded servomechanism whose linear portion is 
described by a complete second-order equation tends to an autooscillatory mode of operation for any parameter val- 


ues satisfying the follow relationships: 


h=4, 





wtg]ch reo 
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Appendix I 





Investigation of the Function v = v(u,y) 
for the Case When h < 1 








The function v = v(u,), for the case when h < 1, is 
defined by expression (14), We determine the signs of the 
derivatives dv/ du, and d*y/ du? by (15). The derivative 
dv/ du, is greater than zero when t, varies in the interval 
(0, ty.) and is less than zero when t, varies in the interval 
(tys, ™ / wy) since 


<0 ininterval (0, qs), 


hi (wf 


>0 in interval (tig */*1). 
fa (ts) <0 


We now determine the sign of the second derivative, 
The sine of the angle w,t; will be greater than zero when 
t, varies in the interval from 0 to ™/w;, The expression 
k-2h+ 1/k is greater than zero, since h < 1, We investi- 
gate the sign of the function f,(t,) [5]. The derivative 
fj (ty) equals 


in interval(0, */w). 


y(t) =h { 2 [« +p) — =| C08 wyty -}- 


A\3 A\? 
H(t SJ me ( tn eS) om 
Since, by the conditions of the problem| Pet > | <i 


then, for cos wt; < 0, the function fj(t,) > 0, and for 
2 
COS Wyt, < 0, the expression [« +p,— 5 | COs w,t; > 0 


and in the undetermined inequality 


2 («a +p,)*— Fos wt V 


(14+m—$) M+(140,44) om 


both members may be.squared without changing the sign 
of the inequality 


sft ug] cost v (1+ 1,4) oh + 


+2[a+uye—F] +(t+m. +5) ™ 
or 


23 
0<4[a+uy'—F] sintoyt, 


+[(14e—F) + (14H4f) em], 


This is, the sign of indetermined inequality V is 
the sign of definite inequality <. 

Thus, fj(t,) is greater than zero for all t;. In the 
‘interval (0, «/ wy) there exists, for hcloseto unity and 


Jet z]<1 such a ty = ty, that f{ty;) = 0 and 


<0 in interval ( 1), 
fo (ts) >0 = in interval (11» */w1)- 


for h close to zero and |e + +| <i, the function 
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f,(t,) is less than zero in the interval (0,7 /w,). Con- 
sequently, the second derivative d*v/ du? can take the 
following signs: 





A 
a { >0at(O,x/w) for hy ~0 |, +> I< t, 
du) >Oat (0, ty) A 
1 for h, ~1 L> <a 
< 0 at (t:,, =/w) Pex D I< 1. 





For t; > 0, the function v > a and u, ~ o, but the deriva- 
tive dv/ du, tends to unity, For t,-* 0, the function v = 
= v(u,) approximates to the iasymptote 


A 


p= 1 K(k — 2h + 1k)’ 


For ty > 1/ wy, the functions v  o and u,*>- o, 
but the derivative tends to -e -™h/w,. The function 
Vv = v(uy), for ty->™ /wy, tends to the asymptote 


(2h —k) (4 +e nh /wy 





on —nh/w ee ae 
a iat k—2h-+- 1k - 
A’ ; A 
— th / wy, ' 
P (uv +5 )e thy 
4 k (k — 2h + 1/k) 


Appendix II 





Investigation of the Function v = v(u,) 
for the Case When h?> 1 








The function v = v(u,) and its derivatives dv/du, 
and d*v/ du? are defined by expressions which may be ob- 
tained from (14) and (15) by replacing wy, sin w,t, and 
Cos wyt, by w, sinh wt, and cosh wt,, respectively, 

We determine the sign of the derivative dv/ du. 
The function f,(ty) is less than zero for 0< t, < ty, and 


greater than zero for ty; < t, < oo, since f,(0) < 0 and, 
for ty = tys, the function f,(t,,) = 0 and f}(ty;) > 0. The 
function f,(t,) is less than zero in the interval 0 < ty < o, 
since f,(0) and f7(ty) = -(1 + y «7 A/ 2)ehty cosh wt, < 0, 
Consequently, the derivative dv/du, is greater than zero 
in the interval 0 < t,; < ty, and less than zero in the inter- 
val ty, < ty © @, 

We determine the sign of the derivative d’v/ du’, 
For t; > .0, the function sinh wt, is greater than zero, The 
function fs(t;), fortyequal to zero, is negative, For t, 
equal to t;,, the function f,(t;), equal to 


[ A? 
(— ch*wt, , + 1) | (i+ p,)?— r 


” + , , 
wehwt, +hshw t,, 


is also less than zero 





and, for ty=oo, equals infinity, Moreover, the derivative 
ff(t,), for t; = 0, equal hw and is positive, Further, 
the derivative f,*(ty) equals zero for ty = ty, and t; = ty, 
whereby 


and 


{ 


A\ 
\t+u+a)e 


—wi ht 
i“. \ 1 ., 


where ty < ty3. Since the function fyt,) for ty = ty, equals 
A? ;' 
[is +, — =| (kh +- w) ee [4 — ehh] and is less 


than zero, then f;(t,;) changes sign only once for some 
ts > tys. In the case when h > 1, the expression 


( k—2h +z) is positive for 0< k< h - dh?-1 or 
k> h+ ¢h*=1, and is negative for 
h—VR—-1<k<h4+Vht~ 1. 


It thus follows from the analysis made that the derivative d*v/ duj_ may assume the following signs: 


— ence or 
h—w<k<h-++w 
O<k<oh—w or 

<0 for alia ea 


For t; ~ 0, the functions v oo and u; ~ o, for 
0<k< h-wor k> h+w, and v*-omanduy*~-o@ 
for h-w < k< h+w, The derivative dv/ du, tends to 
unity for ty->0. Fort; 0, the function v = v(u,) ap- 
proximates to the asymptote v = uy~ A/ k(k-2h + 1/k), 
For t,* oo, the function u, tends to @ for h-w< k<h+w 
and to - .w for 0 < k < h-w or k > h + w, and the func- 
tion 


Jom) + @—w) +(4 —F) = 


k —- 2h + 1)F ° 


A 
afte 





v= 
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k>h+w for 0< <b, 


big < ly < tis ’ 
for tis << 43 < 09; 


k>h-+w for tyr<i<o, 


for O< ty <big bs <i4< hs. 


For t; > oo, the function v = v(u;) approximates to 
the asymptote 


A F A\ 1 
>) (h+w) +(h—w) +(e — ZF 


k—2h-+1)k 


—k—(»,— 
v= 
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Brief descriptions are given of the methods of constructing and correcting signal systems which were de - 
veloped by M., A. Gavrilov [1, 2, 3], R. R. Varshamov [4, 5], and R. W. Hamming [6], The structural formulas 
of the corrector scheme for each of the correction methods described are found. Corrector schemes are contrived. 


One of the methods of providing for correct trans- 
mission of discrete signals entails the use of a corrector, 
which may be connected anywhere in the transmission 
channel, 

We shall assume that the correctors are built of di- 
odes and relay-contact elements, 

Every discrete signal can be presented in the form of 
sequences of zeroes and ones, The number of binary digits 
(bits) in a signal is called its length, There are 2" dif- 
ferent signals of length n, In this set, one defines the 
concept of the distance between any pair of codes a and 
B as the number of noncoincident bits in them, i.e., 


n 
p(a,B)= Dd} o%1 — Bil, 
i=1 
where & = (4, %,..., %,) and B = (B,, By,..., 8p). 

It is well known [6, 3] that, in a system capable of 
correcting d errors (i.e., d simultaneous inverted bits), 
the distance between any two codes in the system must 
be not less than 2d + 1, 

Various methods of constructing such signal systems 
and implementing error correction are known [4-7 ].Here, 
we shall consider only three of them. 

Gavrilov’s method of signal correction, Let a set N 
of working signals be so constructed that any two elements 
of N are separated by a distance not less than D = 2d + 1, 
There is established such a correspondence between the 
signals in N and the set M of all 2” signals that each sig- 
nal « in N corresponds to the same signal « in M and 
that all signals‘ a’¢M, whose distance from a¢M 
is not greater than d correspond to the same signal 
ath. Then, instead of correcting the signal a", one 
can require that each a" be executed as the signal a 
which corresponds to it 

This is a convenient method in that itcan be realized 
for correcting any number of errors d that is desired, 

Hamming's method for signal construction and cor- 
rection, This method is applicable for the correction of 
single errors only, The construction of the set N {s carried 
out according to definite rules, and signal correction (i.e., 
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error detection) reduces to the verification that these 
rules are satisfied, 

Let it be required to construct a set N of working 
signals which contains not more than k signals, There 
is then a number 2 which satisfies the inequality kx 2. 
Then, from the relationship 

J —t—1>1>2°*—(¢ — 1)—1 
t is determined, and the signal length is found to be 
n=1+t, t bits of the signal (specifically, reading from 
left to right, the first, second fourth, eight, sixteenth, 
etc., i.e., the bits in positions 2°=1,2=2,27=4,..., 
2t-1) are called the check bits, while the remaining 1 
bits are the information bits, The information bit posi- 
tions are filled by all the possible combinations of zeroes 
and ones, The check bit positions are filled in accordance 
with the following rules, 

In bit position 2°= 1 there is placed a bit which 
will give a sum of zero modulo two to the bits in posi- 
tions 2q + 1(q=0,1,...). 

In bit position 2! = 2 there is placed a bit which will 
give a sum of zero modulo two to the bits in positions 
4q + 2 and 4q + 3(q=0,1,...). 

In bit position 2? = 4 there is placed a bit which 
will give a sum of zero modulo two to the bits in positions 
8q + 4, 8q + 5, 8q + 6, and 8q + 7(q=0,1. ...). 

In general, in bit position 2! there is placed a bit 
which will give a sum of zero modulo two to the bits in 
positions 2i +1q + 2i, of tig 4 of+1,..., 2b t+ 
+ 2f+1-1(q=0,1,...). 

We now give an example of the construction of the 
set N by Hamming's method, Let k = 4, then 7 = 2,t = 3, 
and n= 5, The signals of set N are given in Table 1. 

It is easily verified that the distance between any 
two signals is not less than 3, 

Error -checking is executed by verifying the rules for 
the construction of set N, All three verifications are made; 
if the results of the verifications are denoted by the letters 
—, © and n, respectively, then the bit position in error in 
the signal is given by the binarynumber € ¢n (if one takes 
1=100,..0,2=010. ..0,4= 00100... 0), For example, 








-~ ow ® «J 








TABLE 1 


























Bit position | 4 2 3 4 5 

Position weight| P BP ¥F > 2? 
ay | a eee 
dy $5? Qe Bly < paleq 
as | “ee eer at ia 
a | yo Eee FG 








if the signal sent were a@,=(11100) and the signal re- 
ceived were @," = (11101), then € = 1, © =0 andy =1, 
The bit position of the error is given by the number & &n, 
equal, in this case, to 101 = 5, 

Varshatov's method for signal construction and cor- 
rection, The set N of working signals is so constructed 
that its elements form an additive group with respect to 
the operation of bit-by-bit addition .modulo two, i.e., 
with respect to the operation defined by the relationship 








a+b=(a,+ by, da + be, oe. , Gn + Un), 


where a = (ay, ag,- ++,» 4p), b= (by, byg,..., bp) and 
the sign "+" is to be understood as addition modulo two, 

Varshamov [4, 5] developed a methodology for con- 
structing such sets N for any given code distance D, We 
limit ourselves here to just a brief description of this 
methodology. 

Let it be required to construct a set N of signals con- 
taining not less than k elements, We determine the least 
integer / satisfying the relationship k = 2! . We then find 
the least integer t satisfying the relationship 


14+ Cart Chat...+ Cra <2. 

Then, n=1 +t, 

The first 2 bits in the signal are the information bits, 
the remaining t bits being the check bits, 

The construction is executed as follows, One first 
constructs a unit matrix of 1 rows and 2 columns, to which 
one appends a second matrix of ! rows and t columns, 
With this, the rows of this latter matrix must satisfy the 
following conditions; the number of ones in each row is 
not less than D-1; the number of ones in eachsum (modulo 
two) of any two rows is not less than D-2 and, in general, 
the number of ones in each sum (modulo two) of any i 
tows is not less than D-i, where i assumes the values i= 
°4,9,.. 0". i 

The set of / signals thus constructed is called a basis. 
The entire set of working signals is comprised of the set 
of basis signals plus all the possible sums (modulo two) of 
any 2,3,..., / basis signals (including the sum of any 
signal with itself, i.e., the null signal), 

Example, Let k = 4 and D = 3, so thatd = 2, We now 
determine t. For D= 3, we get 


Ff 8 {> 


From this, t = 3 and, consequently, n = 5. We con- 
struct the unit matrix and append to it a matrix of two 
rows and three columns, in each row of which there is not 
less than two ones, the sum of these two rows containing 


not less than one *1* 
10 O14 
O1 144 


The set N is made up of the two basis signals, their 
sum (modulo two) and the nuli signal, It is given in 
Table 2 (and coincides with the signals given in Table 1), 














TABLE 2 
ay | 00 000 
a2 | 40 O11 
as | 11 100 
a | O1 411 











R. R. Varshamov developed [5] a method for correc - 
ting d errors if the system of working signals has a code 
distance of D = 2d +1. We limit ourselves here to a de- 
scription of the method for correcting unitary errors, L.e., 
the case when D = 3, 

Let there be received some signal a, 

If follows, from the method of constructing the set 
N, that the information bits of all elements of N exhaust 
all the 2 different sequences of binary numbers, There - 
fore, for each received signal there can be found, among 
the elements of set N, one element whose information 
bits coincide with the information bits of signala, We 
denote this signal by8. We compare their check bits, 

It can happen that their check bits coincide exactly, |.e., 
a = 6, and this means that the transmitted signal was 
received without distortion (error), 

As was shown by the author of the method, if a unitary 
error changes one of the check bits of the transmitted sig- 
nal, then the check bits of signal a will differ by one bit 
from the corresponding bits of signal 8, where the bit 
position of the error is the bit position of the noncoincid- 
ing bit, If the unitary error changed one of the informa - 
tion bits of the transmitted signal, then the check bits of 
signal o will differ from the same bits of signal 8 by 
more than one bit, i.e.,a+B6=y. 

With this, the information bits of y will all be zero, 
and there will be more than one “one” among the check 
bits. To determine the bit position of the error, one seeks, 
among the signals of set N, a signal with the same check 
bits as y (as was shown in [6], if the transmitted signal 
underwent no more than one distortion, then such a sig- 
nal must of necessity be found), We denote itby 6, y 
and 6 differ by just one bit. The position of the latter 
gives the position of the error in signal a, 

We give an example. Let the signal a, = (10011) be 
transmitted, We assume that the signal a = (10010) Is re- 
ceived, i.e,, that an error distorted the last check bit, 
From the table we find the unique signal a, with the same 
information bits as ina, Their check bits differ in one 
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place (bit position), The position of the differing bits is 
the fifth one; consequently, the error distorted the fifth 
bit, which was in fact the case, 

We now assume that the same signal a, was trans- 
mitted, but that the error distorted its second bit, i.e., 

a = (11011), From Table 2 we find the unique signal 

9, = (11100) with the same information bits as a, Their 
check bits differ in three positions, We then compute the 
sum a + B = y = (00111), 

In Table 2 there is a signal 6 whose check bits would 
be the sequence 111, Obviously, 5 = a,, By comparing 
the information bits of y and 5, we convince ourselves 
that the second bits do not coincide, i.e,., that the error 
distorted the second bit, which was in fact the case, 

With some modifications, this method can be extended 
to correct any number of simultaneous errors, 

Synthesis of the correctors, We now determine the 
structure of the corrector for each method of correction 
implementation, 

If each bit of a signal is simulated by the contact of 
some relay, then the signal itself will be simulated by the 
constituent units, 

The corrector schemes are relay devices at whose 
inputs are applied distorted or correct signals from the 
constructed system of working signals and at whose out- 
puts appear only correct signals (the distorted signals will 





A, — Xi Lolgr4rs; 


al —— £1 X%qX%_gX%4X5,; 


alt — ZyXoXygX4rX5; alt — LyloXgX4qX5; all 
a — 24 L2Fgrqr5; _ —2XyX2ryrqr5; a: 
ai’ — £1 LoXgX4rs5; as’ — 14 %22%gX4rs5; a, 


Vv a Vv - 
A, — 2 LeLyrqXs, A, — TyXoLgrgrs,; 


We shall denote the corrector’s outputs by Y, (the 
state of the i th output coincides with the undistorted 
state of the { th input). 

Output Y, must be high (must carry voltage) if the 
second (,) or third (%) signal was transmitted or if any 
of the signals at unit distance from a, or & 3 was applied 
to the corrector input, 

We denote by c, the sum of the constituents corre - 
sponding to a, and e mf aol adh a, , and, by 
Cs, the sum of the constituents corresponding to 


I ll Ill IV 
43, Aes Ox : a, , a, ; and ay. 


Cy == Ty LoF3%g%y + TyTePgT qs -- LyLglgh Ts -| 
4+ LyhghyLyly tT yTeTyXq%y + LyL2T qT ys, + 
+ ¢ = yF oF 37 4s -L 409% 32475 + TF oXgFg75 + 
Hf 2X ot gheFy + XyF 2X pEyTs34 TyFoT37 475, (1) 
which gives us 


Fry, = C2, + C3. (2) 
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a! — LyLolgi4Ts, 


* — Ly LoXgXgX;; a. 


have been corrected by the corrector), It is very import- 
ant in this connection that the number of (simultaneous) 
errors in the signals be matched by the code distances, 

It is obvious that such devices must be nonsequentia], 
With identical systems of working signals, identical sets 
of constituents are applied to the correctors* input, and 
the circuits, in which one and the same set of elements 
operate, can be made identical, Therefore, the correc- 
tors’ schemes can also be made identical, despite the 
fact that the sequences of operations to be executed by 
them are different, Obviously, this assertion is valid for 
any n and d, 

We illustrate its validity by the example wherein 
n=5andd=1, 

Corrector for the Gavrilov method, We write the 
states of relays X;, Xp, X3, X4, and X;, corresponding to 
the signals of Table 2: 


2y— Ty Teng syst 
Ao — Ty VoGgMyls, 
A3g3— 1 LoIyrgX;, 


ae— ryX2 a%gXs. 


We write all the signals at unit distance from ay, 
Og, Xs, and a,, denoting them, respectively, by 


. = Se ae ow 8 ED a 
yy Hy Hy Hy 4 Hy y Ay, My yee rMy s 


a! —— ZyXqXgXqIs; 
sl aisle II a oe 
— LyTq%_gX4rs; a, — TyX2qXgXqr5,; 


— LyLoXgX 4X; gill — LyLeql3gX4rs, 


lV 


— £y%_X3XqX5; Ay — LyLel3l475; 


acs LyLoLgi4r;. 


By introducing the analogous notation for c, 


Cy = LjLoT a7 475 7 &yLo%gT 47s ~< 


T° Fy LolgT gly + Ly Xorg% 4s, +- 


~~ VyLo%gVqFy +- Tj{XoXg%qX, (1") 


and by reasoning the same as in the derivation of Fry,), 
we get . 
Fiyy = C3 + Cy 
Fry, = Co + C1, 

’ 
Foy) = C2 + Ca, (2") 
= Co L C4. 


Corrector for the Hamming method, The corrector 
must execute three verifications, Let the result of each 
of these be recorded by the auxiliary relays Z;, Z2, and 
Z3. Then 


ay ioe 
Fizy=%1E6123¢ Tp = 41%375 


+ 117375 5- LyFgXg + 11 T3Xpy 





or 


ls 
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ch 
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F (z,) = 12 @ Lz = LqXg + Log, 


F(2,) = 14 D Ly = TeXy + Xp (3) 


The circuit comprised of 212, Zs indicates an ab- 
sence of errors, the circuit 2,2 Zz, indicates that an er- 
ror occurred in the second bit, Z,z.z, indicated that an 
error occurred in the third bit, etc. On the basis of such 
arguments, we can arrive at the functional schematic 


shown on Fig, 1. Such a scheme was developed in prac- 
tice and was described in [8]* but, due to the presence 

of feedback in it (from relay Z; to relay X;), its opera- 
tion is unstable, To eliminate the instability of circuit 
operation, one can transform it in the following way: at 
the circuit's outputs we supply the states of the input ele - 
ments and the corresponding corrections to them, obtained 
by means of the auxiliary relay. Then, the action func- 
tions for the input circuits are written as follows: 


Fay=% (22925 T Z\2e2g +> 212925 + ZiZs%3 T 232923) 7 Zy2:2e2s, 
Fy) = 2p (212225 + 212223 + 24.2925 + 24Z2%y + 23292q) ++ Lo2rZarq, 
Fey) = 2p (212223 + 212023 + 212923 + 242923 + 212%) + Tp2sZo%s, (4) 
Fay = % (212925 + 212025 + 22925 + Z12e%y + 2y29%3) + Tehilels, 
Fey, = %p (212925 + 212025 + 2202 + 212925 + Zyh9%q) + ZotrZars.- 
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Fig, 1, Fig. 2. 


The functional schematic of such a corrector is shown 
on Fig, 2. 

If the corrector is constructed in accordance with the 
structural formulas (3) and (4), then the corrector will con 
tain three auxiliary relays, in addition to the five basic 
ones, and will be a two-cycle circuit, 

Since we established that, in principle, the corrector 
is a nonsequential device, we carry out a number of trans- 
formations of structural formulas (4), specifically: we re- 
place the letter z; in them by the corresponding circuit 
paths of contacts x; from (3), As a result of such substitu- 
tions, using the notation of (1), (1"), we obtain structural 
formulas (2), (2"). 

Corrector for the Varshamov method, ‘The noncoin- 
ciding check bits in the received signal a and the signal 
from N with the same information bits as in « will be 
simulated by relays Z;, Za, and Zz, Relay Z;, consequent 
ly, must be switched on if, when signal « is applied, one 
of the following four conditions holds: relays X; and X, 
do not operate and relay X3 does operate; relays X; and 
X3 do operate and relay X, does not; relays X, and X, 





do operate and relay Xs does not; relays X; and Xgdo not 
operate and relay X, does, I.e., 


Fiz) = 2, Zoly + TyTo2y + 4Xo%5 + 142 z. (5) 


By analogous reasoning, we can obtain the action func- 
tions for relays Z, and Zs: 


Fiz, = ZZ, + yXaXq + 242% + TsFa%,, 


Fiz) = XY LnXg + LyLo%5 + Ty X_%y + Ty X_%5. (5") 


If only some one relay Z is activated, then the error 
occurred in the check bits where, if relay Z, ts activated, 
then the first check bit should be corrected, i.e., Ys must 
have voltage on it if relay X, is switched off, and con- 
versely, and all the other outputs must be under voltage if 
the corresponding X relays are switched in, If only relay 
Z, is activated, then the output Y, is in the privileged 
position and , if the single relay Z, is activated, then this 
position will be held only by the output Y;. 

It is easily remarked that an error in the first informa - 
tion bit of any signal will lead to a noncorrespondence in 
the second and third check bits, i.e., relays Z, and Z, must 
be activated simultaneously, It is also easily verified that 
an error in the second information bit of any signal leads 
toa noncorrespondence of all three check bits, i.e., re- 
lays Zy, Z,, and Z, must be activated simultaneously, And, 
finally, if there are no errorsin the signal, then there will 
be no noncorrespondences in the check bits, i.e., all the 
Z relays will be inactivated, 


* In the circuit drawings certain annoying musprints 


occurred; in circuit (a) x; should appear instead of x;, 
and x; should appear instead of X;. 
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All that has been said so far can be written in the form of the following table: 

















ff 242923 the Yi:=7 Y3= 22 Y3= Zs VYo= 2% Y;= 2s, 
2 2at%s Y, = 2 Y2 =T7? Ys = 23 Yo= 2%, Y;= = i) 

If ZyZoz3 then| Y, — x Yo= 2 3= 2s Ye=2%% Ys= Zs 
If 212923 then) Y,; =z 22 Y3= 23 Ye=x% Ys=2s5 
If 212923 the Y;=7 Ys= 2, Y3= 23 Vox 2% Ys= 25 
If 212929 then) Y,;= 2, Yeo= 2 Ys= 23 Yo=% Ys= 2s 

















One can conclude from the table that 


Fy,) = % (212923 + 212025 + 212925 + 242925 + 242025) + 11242025, 

Fy.) = Ly (2,2q2y + 24Z92y + ZyZ%q + 242925 + 242925) + Zy2pZ025, 

Fry, = 25 (212925 + 212923 + 22q2q + 242923 + 229%) + %g2i2oZq, (8) 
Fry.) = %q (24292 + 212925 + 212923 + 212925 + 242925) + 24242975, 


F (Y,) = 7% (25292 + 212923 + 242925 + 232923 + 242923) + p21 2923. 


The functional schematic of such a corrector coin 
cides with the one given on Fig, 2. 

By eliminating contacts of the Z; relays on the basis 
of (5) and (5°), and by using the notation of (1) and (1°), 
we alsoarrive at (2) and (2"), 

The synthesis of the corrector circuit can be executed 
in the following way. We construct the (1,3)-terminal 
network which realizes the conductances Cg, C3, and Cy, 
Then, by means of diode networks, we realize the func - 
tion of Fry,» F(Y;), Fcy,)+ F irs 2 and F(Y;)- 

The corresponding circ shown on Fig, 3. 

In the general case, to synthesis a scheme for any 
n, it is convenient to construct a (1,p )-terminal network T 
with a diode network instead of a (1,n) -terminal net- 
work if, in the system of working signals, there are many 
in which the sign *1° arises in the same bit positions, 

With this, the number of diodes R does not exceed O: 





























Q = 2[2°—(p+ 1)I, (7) 
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whereby R = Q if one of the n functions Fry,), F(Yg)+ ++ + 
F(Yn) for example, Fry))s is the sum of the aforemen- 
tioned p functions, ice., 


Fy) = Ci teat... +c, (8) 
and all the remaining functions are all the possible log- 
ical consequences of applying (8). (Obviously, in this 
case, n-1 = 2P-2),+ 

To convince ourselves of this, we use the method of 
mathematical induction, first noting that, to construct 
each new output from two existing ones, two diodes are 
required (Fig. 4). 


—_ 
_ 


Fig. 4. 


We first consider the case n + 1 = 2P, We convince 
ourselves of the correctness of formula (7) for p = 2(n=3), 
In this case 


Fyy=4+¢, Fwy=ty Fy =¢e 


and 
Q = 2(2’—(2 + 1)] = 2(4—3) =2, 


which is in accord with the actual facts (Fig. 4). 
We now assume that (7) is valid for some p > 2, and 
we prove its validity for p + 1, 





T Here, p is the number of different Boolean functions cy, 
Cy,..., C, entering into the expressions Fry, (i = 1, 2, 
« «+» 0), each of which is the sum of the constituent units 
corresponding to some working signal and all of its dis- 
torted variants, 

+We exclude from the number of all the consequences 
the two trivial ones; 0, and the function Fry ) itself. 





of 


and 





Let there be given p + 1 outputs, and let it be nec- 
essary to obtain from them all possible outputs by inter- 
connecting the primordial p + 1. We assume that, from 
p outputs,we have already obtained all their possible 
outputs (their number will be Ch +e e@e...4G@) 
and that we have used 2[2P-(p +1)] Fiodes for this, We 
shall connect the (p+ 1)th output with each of these out- 
puts, using two diode elements for each of these connec - 
tions. With this, we shall use 


25Cch + Cc? of C3 +... CB] = 2? (2? — 1) 


diodes, 
Consequently, we shall use a total of 


2 (2? — 1) + 2[2? —(p + 1)] =2 (2°41 [(p + 1) + 4) 


diodes, q.e.d. 

If n + 1 < 2P, then the number of connections and 
the number of diode elements used will be smaller. 

It should be mentioned that the assumption that the 
correctors be realized by relay contact circuits is not an 
essential one, In [9] there was . presented a method of 
synthesizing correctors from contactless elements by 
structural formulas of the form of (2) and (2°). 
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A description is given of a 10-channel automatic electronic-relay optimizer. The operational part of 
the optimizer consists of electronic operational amplifiers, while the controlling portion is comprised of relays 


and stepping switches. 


The optimizer is designed to operate with electronic analog computers and is used to obtain solutions to 
complicated variational problems, and problems of finding extrema. An example of problem solution by means 


of the optimizer is given. 


Certain problems in finding extrema cannot be solved 
analytically, either because analytic expressions of the 
functions are not known, or because these expressions are 
too complicated, This is particularly the case for the prob- 
lems related to the finding of extrema of functions of 
many variables, In actual systems, or in systems simulated 
on analog computers, one can, by giving increments to 
the individual variables, easily determine partial deriva - 
tives and, consequently, gradients of a function. By mov- 
ing either in the direction of the gradient or in a direct- 
ly opposite direction, one can find a maximum or a mini- 
mum of a function, The processes of automatically deter- 
mining partial derivatives and of moving in the direction 
of the gradient can be implemented by means of the so- 
called automatic optimization systems. If, from the con- 
ditions of the problem, certain boundaries are imposed on 
the variables then, as was shown in [1], the determination 
of an extremum lying on a boundary is possible in this 
case, 

The theoretical considerations bearing on the auto- 
matic search process were presented in [1, 2]. A complete 
electronic variant of an optimizer was described in [3]. 

Below there is given a description of an electronic - 
relay instrument which executes automatic search for ex - 
trema, For brevity, we shall henceforth call this instru- 
ment an optimizer, Participants in the design and adjust - 
ment of the optimizer were I. N. Bocharov, A. V. Kalinina, 
R, I, Stakhovskii, and the author of the present paper. 

The optimizer permits one to find the extremum of a 
function of many variables (up to 10) with the presence of 
limitations (boundaries), 

The function whose extremum is to be sought, as well 
as the function which {s bounded, are applied to the op- 
timizer in the form of de voltages whose range of varia- 
tion is +100 v, The optimizer can operate with objects 
which continuously output the value of the function whose 
extremum is being sought, as well as with objects which 
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output the function value only at certain moments of time 
(for example, after solution on an electronic analog com- 
puter). 

The optimizer can be divided into two parts: the 
operational part and the controlling part, The basic com- 
puting portions of the automatic optimizer are the elec- 
tronic analog blocks, The output voltages x,, Xg,.~ « »Xjg 
are formed on de integrators; the range of variation of 
the output voltages is +100 v. The size of the trial step, 
applied for the determination of the partial derivatives, 
is regulated within the limits of 0 to 15 v. The choices 
just cited are related to the circumstance that the optim- 
izer is basically designed to operate with electronic ana- 
log computers and, with such choices for the computing 
(operational) part, matching of the optimizer with ana- 
log computers is executed with extreme simplicity, The 
accuracy of the operations effected by the computer and 
the optimizer is of the same order of magnitude in all 
cases, The accuracy with which the extremum is deter- 
mined depends on the character of the problem to be 
solved, For a test object, which was the sum of the moduli 


10 

of the ten variables Q =k be jay, the maximum rela- 
i=1 

tive error in determining Q.,,, was 1%, 

The optimizer's controlling portion analyzes the data 
applied to the optimizer and analyzes the state of the op- 
timizer itself and, as a function of the information ob- 
tained, carries out some switching or other, Relays and 
stepping switches are the basic elements of the controlling 
portion, Adequate operation of such a controlling block 
is possible for frequencies of the synchronizing pulses up 
to 5 cps, i.e., it requires about 0,2 sec to apply the trial 
increments and to remember (store) the partial derivatives. 
Such a speed of operation is completely sufficient for op- 
eration with ordinary electronic computers which are also 
relay-controlled, 
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Fig. 2, 1 is the block for forming limitations, 
2 is the object, 


An exterior view of the optimizer is shown on Fig, 1. 
The optimizer design incorporates, in one chassis, the con- 
trolling block, the limitations block, the AQ block, the 
test block and 10 channel blocks, The front panel of the 
optimizer is comprised of the front panels of the individ- 
ual blocks, on which are provided the necessary devices 
for manual control and testing of the instrument, The 
dimensions of the optimizer, without the power supply, 
are 650 X 500 X 450 cm, 


Operational Portion 





Figure 2 shows the functional schematic of the optim- 
izer's operational portion, The theory of operation 





amounts to the following: the quantityQ, whose extremum 
is sought, is applied to the AQ block, This block remem- 
bers the function and determines its variation, We assume, 
that, in the AQ block, there has been remembered the val- 
ue of the function for some values of x;, Xg, . . « » Xygi if 
we now give a trial increment Ax;,, to the first variable 
then, at the output of the AQ block we obtain the incre- 
ment of the function AQ, ™ (8.Q/8 x) Axyy, f.e., a quant- 
ity which is proportional to the partial derivative of the 
function with respect to the first variable, Simultaneous- 
ly with the application of the trial increment, contact 
KM, Closes, and storage of the quantities #(9Q/0 xjAxy 
occurs in the memory link of channel 1, 

After this, a trial increment is applied to the second 
variable, etc,, up to the tenth variable, At the end of the 
cycle, there have been stored in the channels* memory 
links the quantities which are proportional to the partial 
derivatives, or components of the function's gradient, Up- 
on simultaneous closure, at time ty, of contacts Kl, Kl,, 

. « » Kp, the integrators receive the increments 


Ary =+ “e Ary kts , 


Azz, = + a Az, kt, ete. 


OX 


The trial increments Ax, and the integration con- 
stants k are identical in all the channels, so that the 
operating variations of the variables (the working steps) 
are proportional to the components of the gradient, After 
a working step, one can repeat the determination of the 
gradient and carry out a new working step along the new 
gradient; in this case, the search is carried out by the 
gradient method, After the working step one can deter- 
mine the sign of the function's change and, if this change 
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Fig. 3. 


occurred in the proper direction, then the working step 
along the previous gradient may be repeated until such 
time as the increment changes sign, i.e., until such time 
as a partial extremum in the direction of the previous 
gradient is reached, This method of search is called the 
method of steepest descent, There are a number of 
switches on the controlling block whose positions corre - 
spond to the various search methods, 

If limitat.ons are imposed on the system's coordinates 
i.e., if some function Hgx;, Xg, . . « » Xyg) Must not ex- 
ceed a given value, then the search for an extremum is 
executed in the following manner. Optimization of the 
function Q proceeds until the function Hy goes beyond 
the limits of the admissible region; with this, there is 
automatically switched to the input of the AQ block, | 
not the function Q whose extremum is sought, but the 
bounded function Hg, In this case, x3, Xg,. . « » X49 are 
so changed that they move out of the region of limita - 
tion, After the region of limitation has been left, the 
process of finding an extremum of function Q is continu- 
ed, If an extremum is found either on the boundary or 
inside the region of limitation, there then occurs an 
oscillation about the boundary and, in the process of these 
oscillations, the optimizer finds the true extremum or a 
conditional extremum on the boundary. 
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We now consider the circuits of the individual blocks 
of the optimizer. Used in the optimizer blocks are stand- 
ard cathode followers (Fig. 3 b) and de operational ampli- 
fiers (Fig. 3a), The amplifiers operate in either the am- 
plification or the integration modes, 

In the AQ block, the storing of the function whose 
extremum is sought and the obtaining of its next incre- 
ment occur on condenser C, (Fig. 4a). When contact 
Rin 1s closed, the condenser charges up to the potential 
Q, while, when the contact is opened and Q varies, the 
following voltage is formed at the output of the cathode 
followers: 


AQ = (2 — Qo) hres 


where k,,, is the cathode follower's transfer coefficient. 
The function increment obtained from the cathode 
follower is amplified and applied to the channels’ mem- 
ory links, There are inverters at the input of the AQblock 
which are switched in to the input Q by toggle switch 
To and to input Hy by toggle switch T,,,. This makes 
it possible to change the signs of the partial derivatives 
@ Q/ 2 xiAx;, and (0H, /2 x;)Ax;,, and thereby to carry 
out either minimization or maximization, 
Figure 4b gives the functional schematic for one 
channel, Storage of the partial derivative occurs on con- 
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denser C, upon operation of storage relay Rm, 1.e., upon 
closure of contact KM. In order that the charge on the 
condenser not be decreased when connected in to the in- 
tegrator’s input, the voltage is taken off from it by means 
of the cathode follower. The cathode follower is based 
on tube 6 N 3P which has low grid current; the change in 
the condenser’s charge due to grid current does not exceed 
0,1 v/ min. To eliminate potential leakage from con- 
densers C, and C; (and also from condenser Cy in the 
AQ block) as well as a fatigue effect, polystyrene MPGT 
condensers of 1 yf capacity are used in the optimizer. 

The output voltages x;,X9, ... » Xy9 are supplied 
from the de amplifier-integrators, The amplifier has 
high gain, uses tube 6N3P in its first stage, and the operat- 
ing mode of its first stage is chosen so as to minimize 
grid current, Integrating capacitor Cy is connected in the 
amplifier’s feedback path. 

Application of the trial step is implemented by means 
of condenser C,, Upon closure of contact RR at the integra - 
tor input, a step voltage Urs is applied via the condenser, 
As is known, the integrator is an amplifier under these con- 
ditions, so that the following increment appears at the in- 
tegrator output, 


Ar, = Ths i . 

Upon opening of contact RR there is a voltage step 
of the opposite sign, and the increment Ax; is tapped off. 

The working changes of the output voltages x,,X», 

. « » Xy9 (the working steps) occur upon operation of step 
relay Rs, i.e., upon closure of contact KI. With this, the 
integrator input, via a resistor of 1 meg, is connected to 
the output of the cathode follower, the voltage on which 
is proportional to the partial derivative with respect to 
these variables, Thus, the output voltages receive incre - 
ments proportional to the partial derivatives, 

The initial conditions on the integrators are estab- 
lished by connecting toggle switch T, to the channel 
block, With this, the capacitance of the feedback path 
is shunted by a resistance and, at the integrator's input, 

a voltage from a divider is applied. 

The circuit for forming the limitation Hg is given in 
Fig. 4c. The circuit includes an amplifier with four in- 
puts, There is a vacuum-tube diode at each input. Thus, 
the circuit gives the sum of only the positive voltages 
Hj which is just the bounded function Hg. The function 
to be limited is applied to the limitation indication cir- 
cuit in the indicator block. 

The circuit for indicating limitiations (similar to the 
circuit used for indicating partial minima and indicating 
the final solution) is shown on Fig, 3c, The circuit works 
as follows, In the differential stage’s anode circuit there 
are connected the windings of the polarized relay RP-4, 
The windings are connected toward each other, On one 
of the grids, Ge, the constant voltage Us, is applied, while 
on the second grid, G;, the voltage U to be compared is 
applied, If U > Up, then the current in winding L, is 
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Fig. 4. 


greater than that in winding L,, and the polarized relay 

is found in one of its positions, For U < Us, the relay is 
transferred to its other position. Such a circuit allows 
one to fix the transition of a voltage through any threshold 
level with an accuracy of about 0,5 v; the hysteresis of 
the circuit is less than 0,1 v. 


Controlling Block 


The controlling block is comprised of a stepping 
switches, relays and switches of various types of opera - 
tion (Fig. 5). 

The first stepping switch SSI executes successive 
switching in of the memory relays of the channels and, 
by the Gauss-Seidel method, also executes the individual 
switching of the channels’step relays, Simultaneously 
with the switching of the channels* memory relays, the 
trial increments for the variables are supplied, Thus, as 
SSI moves along commutator segments 2-11, the gradient 
is determined for the function whose extremum is being 
sought, When on the first commutator segment, SSI ap- 
plies voltage to SSII,causing execution of the operations of 
working-step implementation and determination of the 
change in the function during the working steps, With the 
Gauss-Seidel method, these operations are possible on seg~- 
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Fig. 5, 


ments 2-11, while commutator segment of stepping switch 
SSI is passed over by shunt running, 

Various relays implement the following basic func - 
tions: 

1, Start relay Ry, when the start button is pushed, 
switches in cycle relay Re. 

2, Cycle relay R, operates upon the reception of a 
sync pulse from the synchronization block, and transfers 
SSI or SSII to the following commutator segment, i.e., 
implements transfer from one operation to the following 
one, For the optimization of continuous objects, the re- 
lay is transferred to the autooscillatory mode by toggle 
switch S,. 

3, Relay R,, for shortening the variables cycle,switches 
SSI to shunt running in order to pass rapidly over the com- 
mutator segments which correspond to unused channels, 

4, Step relays Rp and Rpg provide the required de- 
lay to the working step, This delay is determined by con- 
denser C, and resistor R,. 

5. The minimum relay R,, operates when a partial 
(relative) minimum is reached, After the storage opera - 
tion, relay Ry,, by its own contact 2Ry4, switches SSI to 
gradient determination, The operation of gradient deter - 
mination must also be executed when the instrument is 
started and upon input or output of limitations, In the 
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first case, relay is switched on by a normally open 
contact of the start button and, in the second case, by a 
special contact 6R,,, of limitation blocking (holding ) re- 
lay Ryyp- This contact is closed only in the periods when 
the relay is transferring from one state to another. 

6. Holding relay Ry operates after the termination 
of a working step, At this moment there is implemented 
the analysis of the increment of the function whose ex- 
trumum is being sought, contact 3R,_ is closed and the 
minimum relay is switched to the circuit for minimum 


indication, 
Contact 1R), releases relay Ryy,, and contact 2R, 


switches the limitation holding relay to the limitation 
indication circuit, After release of R,, the state of relay 
Rip depends on the state of relay Ry at the moment of 
holding, and remains unchanged until the following ope- 
ration of relay R,. 


For the optimization of functions obtained after a 
solution on an analog computer, it is necessary that the 
holding relay operate at the moment when there is a val- 
ue of the function to be optimized available at the com- 
puter, For this, 2-holding is used, 

For an extremum search using the Gauss-Seidel meth- 
od, the memory relays in the channels are connected-in 
via the commutator segments of SSI and contact 3Rws 





ay 





of the working step relay. To switch this relay in, there 

is auxilfary relay R,. Relay Ry operates simultaneously 
with SSI. Thus, when SSI transfers to the following seg- 
ment, i.e., when the optimization cycle for the previous 
variable terminates, relay Rx is found in the operated state, 
Since contact 4R, is closed, then simultaneously with the 
operation of Rp working step relay R,, operates and, by 

its contact 3R,,., it connects in the memory relay in the 
corresponding channel; the operation of partial derivative 
determination then ensues, Contact 5R,, cuts out the 
channel working step relay, Simultaneously relay Ry, 

by contact 4R,., cuts out relay R,. At the time of the 
following working step, contact 4R, will be open, relay 
R.,, will not operate, so that a working step will be ex- 
ecuted, etc, With the Gauss-Seidel method, the channels’ 
working step relays are opened by switch S3, and the work- 
ing step is executed for one of the channels. 

The switches have the following basic functions: T,; 
applies voltage to the circuit; the closed position of 
switch S, corresponds to search by the gradient method, 
i.e., after each working step the determination of a new 
gradient is made; S, switches of optimizer from opera - 
tion by the gradient method to operation by the Gauss- 
Seidel method; S, switches in 2-holding (secondary hold- 
ing) for discrete objects; S; connects the cycle relay to 
a source of external pulses; S, reduces the number of 
variables for optimization. 


Solution of a Problem by Means of the 
Automatic Optimizer 








In addition to the solution of problems without limita- 
tions, such as the automatic approximation of a given 
curve [3], the optimizer permits the solution of complica- 
ted problems with limitations, 

The problem arose of finding the optimal control pro- 
gram for a motor-generator system such that, during a 
given time T, the motor would make the maximum num- 
ber of rotations from the moment of starting to the mo- 
ment of stopping, while observing the limitations imposed 
on the system's internal coordinates, 

As the first step, it was assumed that the motor's ex- 
citation current remains constant, An approximate linear- 
ized system of equations for the object has the form: 


dF. 
Tp + Eg= Wi, 


E,— cQh = L Ra 9 < = FQ. 

Here, U, is the voltage applied to the generator's 
excitation winding, Ey is the generator’s emf, I, is the 
motor’s armature current, I.,, is the motor's excitation 
current, Q is the angular velocity of the motor shaft, X 
is the angle of rotation of the output shaft, connected by 
a reducer to the motor shaft, R, is the resistor of the cir- 
cuit of the motor armature and the generator ar™ature, 

The circuit set up on the analog computer to simu- 
late the given differential equations is shown on Fig, 6a. 


a 


u,(t) 














Fig. 6, 


The program for varying the generator's excitation 
voltage was supplied from a special synchronizing block. 
The program was divided into ten identical time inter- 
vals and, in each interval, the voltage was chosen by the 
optimizer, 

The following limitations were imposed on the sys~- 
tem: 


|U,|<U, max |, |<Ja max |Q|<Qmax 


In the given problem, it was required that the final 
conditions Qeng = 0 and E, = 0 be met, These conditions 
derive from the posing of the problein, since the complete 
stopping of the motor at the end of the process is required, 
For the meeting of these conditions, the following addi- 
tional limitations were imposed: 


Qend< Nandy end <0, 


After starting of the system from null initial condi- 
tions, the optimizer automatically found the optimal pro- 
gram. In the majority of cases, the maximum X lay at 
the point of intersection of the boundaries of limitation: 


Qend< 0 Es end 59% I; <1, max: 


The optimal processes found by the optimizer are 
shown on Fig. 6, b and c, 


CONCLUSIONS 


The automatic multichannel optimizer is an instru- 
ment allowing one to solve variational problems and 
problems of finding extrema of functions of many vari- 
ables, 

The given optimizer was designed basically to operate 
with electronic analog computers with relay control, How- 
ever, the principles on which the instrument is based can 
be used for creating optimizers to operate with actual ob- 
jects, It is possible to have optimization of actual ob- 
jects carried out in two time scales, iie., rapid finding of 
extrema on an analog computer and a subsequent transfer 
of the results to the actual object. 
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Methods are considered for the presentation of reliability characteristics of resistors and condensers as func - 
tions of ambient air temperature and power dissipation or applied voltage. 


The reliability of constant carbon resistors, among the 
most widely used elements of different systems, depends 


on the physicochemical processes occurring both within the 


conducting layer and on its surface. Therefore, the relia- 


bility of resistor operation is determined, to a large degree, 


by the choice of materials, the design, and the technology 
of manufacture, In addition, the reliability of resistor 
operation also depends on thecorrectness of the selection 
of its operating regimen and on external conditions, 

In designing circuits containing resistors and conden- 
sers, one must know how to determine their reliability un- 
der given conditions of use, which may be characterized 
by the ambient temperature 6x, dissipated power P,, volt- 
age U, and length of serviceable life T. For this end, it 
is convenient to introduce characteristics which give the 
connections between reliability and the parameters which 
define the conditions of use. The consideration of these 
connections comprises the subject matter of the present 
paper, 

1, By considering the breakdowns of resistors to be 
random events, we have that the number of operational 
failures during time interval dt will be the larger, the 
greater is the number N of resistors, i.e., 


—dN =cNdt. 
It follows from this that 


(1) 


InN =—ct+Ink or N=ke-“, 
Since N = Ny = k for t= 0, then 
N = N,e-* “4 


or Ne = p = e~ct” 
Here, N is the number of existing elements (resistors) at 
time t:N/No=p is the reliability of action of the set of 
No elements at time t; c= 1/r is the reliability factor: 
(the relative speed of breakdown), 

Experimental data substantiate the foregoing depend- 
ence of element reliability on time for the fundamental 
portion of the characteristic p = N/ Ng = f(t) (segment AB 
on Fig, 1), Deviations occur only in the initial period, 
when the unconditioned elements go out of order, and af- 





* 
‘init 
Fig. 1 


t 
op. lim 


ter point B when, due to accelerated aging and to other 
irreversible changes, the element breakdown rate begins 
to accelerate. 

Point B determines the length of serviceable life T 
for the given conditions of operation, 

We now consider the dependence of the coefficient 
c on the resistor's heating temperature 6 ,. 

The velocity of change of material during the reac- 
tion process is defined by the equation 

dn 
eden 

where n is the quantity of material at the moment of time 
under consideration and a is the constant velocity of the 
reaction, 


From the equation given it follows that 


oo —adt, 
n 
i.e., 
T 
lhna=— \ adt, 
0 
or, for « = const, 
In> <= —aT, 
n, 


0 
The dependence of the constant reaction speed on 
temperature is determined, according to Arrhenius, by 
the formula 
a= A, exp(— a ’ 
. 273° + 0, 
where A; and B are constants and (273 + @,) degrees is 
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the absolute temperature, 
Consequently, 


=e)? 


n 
In a A, exp (— s+ 0, 


We find from the last expression that 


T=|5 In “*| exp (agra) 


or 


T = To. exp (ae a) -)s 
If 6, $275°C, then 


exp ( 
T 


, ga B 
oo (ors * 8 a) 


B (6, — 9%) 
(273° + 6.) (273° + 0) .) 


B ) 
273° + 6, 





a5 





= exp a 


If 6. < 273°C, then 
B 
T a T co exp (a - a.) = 


B 
= T.. exp = 
273°( 6, \ 
D7 


B B 
= T., eXP 579 EXP (— 738 6)= 
= T, exp (— £6,), 








or 
T == Tyexp([— B (0x — %)], 
oe = 7, exp YP 
8 = saa » T'y>=T, exp (— £6). 

By denoting by s, the amount of material in one re- 
sistor, we get that n = syN and, consequently -sdN/dt = 
= sgNa or —dN/dt = Na, 

By comparing the last equation with Eq. (1), we get 
that a = c, We denote by cy the value of c which corre- 
sponds to @»9; we then find that 





«1, exp ( — ———_—— 
— B a 
ty A, exp [- FTG) 
B (0, — %) \ 





= exp \ oa -+- 6.) (273° -- Oo) | (2) 


The value of the temperature 6, is determined as the 


sum 
),. = Do. + RP x, (3) 


where @ox is the ambient temperature, R, is the thermal 
impedance for heat transfer from the surface to the am- 
bient medium and P, is the power dissipated in the elec- 
trical impedance, 

By substituting from (3) to (2), and by denoting by 
Cy the value of c for the nominal load P,, and the ad- 
missible temperature 69x = 6», we get 
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P 
B\0..—0.,+R,P | —1] | 
- | CN t \Py } | 
— = 6) ° > 972° 
(213° + 8. + RP) (279 F OF RAY) * 
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Fig. 2 


Figure 2 shows the function c ley, for one type of re- 
sistor, This function is also the characteristic c/ cy 

= f(@ 9x, Px) for the resistor. 

The quantities B, R, and Cy, are obtained by process- 
ing experimental data. The value 69 ) corresponds to 
the maximum permissible temperature for the nominal 
load P, = Py. 

It should be kept in mind that the values of B, Re, 
and c,, vary somewhat as functions of the type and the 
nominal resistance, However, for determining resistor 
reliability in the process of designing a scheme, it suf- 
fices to estimate the average value of reliability for a 
given type of resistor, Such averaged values of c,, and B 
are given in the following table, 

If the limiting admissible value for the resistor tem- 
perature is 6, max, then the magnitude of the admissible 
power dissipation equals P = (8x max~ 8 ox )/R,. 

Figure 3 shows the lamnhtbe P. = £(6,). 

Thevalue of @, = @nom = 6 xN is chosen from the 
condition limiting the relative variation of the resistance 
AR/R= f(6,) for ox = 6 ny» and should be a value which 
is achievable technologically or which is determined from 
the given conditions of use. 

The method of choosing the value of P,.., = Py con 
sists of this: that knowing the length of useful life Ty and 
the temperature 6, corresponding to it, we give a value to 
the admissible service life T. Then, based on the fore- 
going, we obtain 


1 B © oS) 
Sat 9 Se Sempre steed Lt Mia Os) : 
eerie + 18 
273° + 06 y. j 
With the condition that the value of c remain invariant, 
the limiting value 65, = 6x max is determined if we 
set P =G6, i.e., 
x 
B 
B ‘ 


a a) 
273° + Oo To 





0, max >= Box = — 273°. 
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Nominal|Admissible | Limiting [Thermal C), 
Type of power /limitin ambient |; 
g impedance (mean 
resistor in W om temperatura, . 8 
perature) 6 (for p = [Re in value), 
®max in°C | By), in °C |*Chw hr 
VS —0.12(ULM—| 0.12 130 180 350 
(), 12) 
S —0.25 0.25 130 95 140 
VS —0.5 0.5 130 85 90 1.0-10-8 |0.13—0.18 
VS —1.0 1.0 130 65 65 
S —2.0 2.0 130 70 45 
S —5.0 5.0 130 35 19 
S — 10.0 10.0 130 50 8 
S per spec. 0,12—10,0 110 40 70 
Lt —0.5 0.5 130 70 50 
— 1.0 1.0 130 _ — 1.7-10-* |0.11—0.12 
LT — 2.0 2.0 130 —- apes 
=A@. The relationship can be expressed approximately 
P, by the formula 
Pp, V1 = 10 (A6)™, 
where 149 Corresponds to the temperature differenceA 6 = 
— = 1°C, 
p.! > Thus, the value of the nominal power is determined 
er from the expression 
0 657 9.0m% yn 6, a * Ph 
ax a, ma max, Ri =(0, sax — 9) Ss Uyo(A0)™ _ - 
Fig. 3 Fig. 4 : 


Conversely, by giving various values: c = cy, C =C3, 
C Cs... , we obtain various magnitudes of 6, maxe 
This makes it possible to define Py for a given @y = const 
(Fig. 4). The graph of Fig, 4 is another form of the same 
characteristic c/cy = f(@9,, Px) which relates the change 
of reliability to the change of temperature of the ambient 
medium and to the change of power dissipation P,. 

The value of the thermalimpedance can be found from 
the formula 





4 

Ki = SF UST1’ 
where yp is the total heat transfer coefficient, S is the 
cooling surface, A is the thermal conductivity of the 
current leads and S» and/ 9 are their cross section and 
length, The expression for the heat transfer coefficient 
is written in the form 
(273° + 6,)* — (273° + 6,,)* 

6. — 6. 





B=u,+5 ’ 

where p, is the heat transfer coefficient due to gas heat 
conduction and convection, It should be kept in mind 
that the value of p, varies with changes in p, _, the pres- 
sure of the air surrounding the resistor, by the law 


Ph 
pA = —, 
ih, = Pin, Pr, 
where pp, is air pressure at sea level, The magnitude of 
}, also dependson the temperature difference 6,~6,, = 








c  — ° 
4 ye +9 mad —ar+e "|as 
* max *N " 

It follows from this expression that the value of the 
nominal power must vary somewhat as a function of the 
air pressure and the value of 0) = 6, max~Oy- It will 
also depend on the mounting, since /9and p walso vary 
with this. 

The foregoing table gave the values of R, and cy, for 
a number of standard resistors, and also gave the permissi- 
ble temperature of the ambient medium, the value of 
which changes for resistors of different power ratings, 

2. In the case of condensers, it is necessary in esti- 
mating the reliability of the finished articles, to take into 
account both the effect of temperature and the effect of 
voltage, 

The useful life of paper condensers, determined by 
the aging of the insulation, as a function of the tempera - 
ture 6, is expressed by the formula 


T'xq = T) exp [— Bn (8x — %)}, (4) 
where 6, is a coefficient which depends on the form of the 
insulation, 

On the other hand, for dielectric aging in the pres- 
ence of an electric field, we have the empirical relation - 
ship 


e= AT = ea - (5) 
where € is the electric field strength, T is the length of 
useful life and A is a constant. 
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By taking into account that U = €d we get that U = 
= AdT™, where d ts the thickness of the dielectric, For 
condensers with paper insulation, n ® 0,2, 

For two different voltages U, and Ux» we have that 


( U, ¥" in 

a pe 
In order to take into account the influence both of 

U, and of 6,, we write the ratio Tyo/ Tx in the- form 


Tyo _ To Tso 


i i 





J 


Then, 
U. 1/n T 
(a) = p exP [— Bn (Ox — %)I, 
i.e., 
U T 
ng = —n In 7 — nBn (6. — 6%). (6) 


If the number of condensers remaining in service for 
Ux = const and @, = const is defined by the relationship 


N ’ 
P= H, = xp (— Cel’) 


or 


N 
In p= Inge = — Gls, 
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then, by giving some value of p, one can easily deter- 
mine the value of length of useful life which corresponds 
to it: 


T= ——, (7) 
k 


By substituting the value of T, in (6), we get 





U - 
mt —nin{ — “OP. | — nn (Ox - No). (8) 
xo “hr 0 


Ur 





@,¢ 


Fig. 5 


Figure 5 gives the given characteristic, defined by 
(6), allowing one to relate the length of useful life T,, to 
the admissible values of the voltage U, and the tempera - 
ture 6,. 
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ON THE DESIGN OF MAGNETIC AMPLIFIERS 
SINGLE-PHASE AND THREE-PHASE SEMICONDUCTOR 
RECTIFIERS ON AN ACTIVE LOAD ANDOPERATING WITH 


A COUNTER-EMF 


A, G, Zdrok 
(Moscow) 


Translated from Avtomatika i Telemekhanika Vol, 21, No. 5, pp. 639-651, May, 1960 
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Methods are presented for the transfiguration of single-phase and three-phase rectifier bridges to equivalent 
loads, and graphicoanalytic methods are considered for determining the characteristics of saturated reactor mag- 


netic amplifiers operating on active loads and counter-emf's, 


Circuits for automatic current stabilizers in charging devices are given, and anexample of determining a 


magnetic amplifier's basic parameters is considered, 


Semiconductor rectifiers, in conjunction with mag- 
netic amplifiers, allow the construction of compact auto- 
matic control devices, Such devices are widely used for 
controlling the speed of a motor's rotation by varying the 
voltage applied to the armature, for charging storage bat- 
teries, etc, [1-4], Setups with magnetic amplifiers and 
semiconductor rectifiers are increasingly superceding elec- 
tronic and ionic tubes with grid control. 

The currents and voltages at the input and output of 
rectifier circuits are different, in the majority of cases, 
The relationships between them vary with changes in the 
circuit's operating mode (sinusoidal) voltage or sinusoidal 
current at the rectifier input, etc.) and also with changes 
in the form of the load (active load, inductive load, coun- 
ter-emf, etc.), It is usually quite difficult to find the re- 
lationships between the input and output currents and volt- 
ages, as well as to design the magnetic amplifiers operat- 
ing in such circuits, 

In this paper we present graphicoanalytic methods for 
determining the characteristics of saturated reactor mag- 
netic amplifiers which operate with single-phase and 
three-phase semiconductor rectifiers, these methods be - 
ing based on the transfiguration of the rectifier scheme to 


a load which is equivalent to it, All other rectifier schemes 


can be dealt with similarly to the schemes to be consider- 
ed here, 


Magnetic Amplifier Operation with Single- 





Phase Rectifier Bridge 





With operation of a magnetic amplifier on an active 
load, supplied with rectified current, the instantaneous, 
average and effective currents at the input and output of 
the rectifier bridge are all equal to one another, as are 
the corresponding values of voltage, These relationships 
are considered for the case of ideal rectifiers, for the ab- 


solute magnitudes of the instantaneous values and the aver- 
age values over a half-period, If the voltage at the recti- 
fier inputs is sinusoidal, then the relationship between the 
voltage's average value and effective value, as well as 
the relationship between the analogous current quantities, 
are constant and equal to 0.9, Therefore, the design of 
such a Circuit engenders no difficulty. 

Significant difficulty is engendered in the design of 
a circuit for rectifier operation on a counter-emf since, 
in this case, the relationships between the mean and ef- 
fective values of current and voltage differ from each 
other, and do not remain constant under variations of the 
circuit’s mode of operation (for example, for sinusoidal 
voltage and sinusoidal current on the rectifier inputs), 

Figure 1 shows the variation of voltages and currents 
for the operation of a single-phase rectifier bridge on a 
counter emf for the modes of sinusoidal voltage (Fig, 1a) 
and sinusoidal current (Fig. 1b) on the rectifier inputs. 
Here, u and i are the instantaneous values of voltage and 
current on the rectifier inputs, ug is the instantaneous 
value of voltage at the rectifier bridge*s output and E is 
the value of the counter-emf. 
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In the mode of sinusoidal voltageon the rectifier in- For a sinusoidal current at the input of a single- 
puts, the angles @, and @, are defined by the following re- phase rectifier bridge, the coefficients 6, and 6; have 


lationships [1): the following values: 
, E 
~* Tn” di: eR ies @) 8, = 0,9 for 7 = 0, (9a) 
where Upp is the amplitude of the voltage applied to the Pas 
rectifier bridge. | B.=1 for W- = 4, (9b) 
We denote by 8; the ratio of the mean value of recti- g —e 
fied current Ig to the effective value of current I at the B,=9 | 0< v- <1}. (9c) 


input of the rectifier bridge, i.e., ; 
Due to the inconstancy of the coefficients 8,, and 


te V4 (2) B; for a variable ratio of E/ Um, and also for sinusoidal 
glad current and sinusoidal voltage at the rectifier inputs, it 
and, by 8,,, the analogous relationship for voltage, i.e., is much more convenient, in designing the circuit, to use 
U4 the replacement circuit with an equivalent active load 
by = T° (3) re instead of the rectifier bridge. 
With a sinusoidal voltage at the input of a single - For operation on a counter emf we have 
phase rectifier bridge, the coefficients 8; and 8, have Us=E+ lyr, (10) 
the following values [1): 
E where r is the active impedance in series with the coun- 
8; = 0.9 8. = 0.9 for =— = 0 (4a) -emf E (i . 
Pi , u U ’ ter-emf E (if the rectifier output is connected to the ar 
: os mature circuit of a motor or to a storage battery, thenr 
8B.=4.44. 8,—1 for —wzi. (4b) is the active impedance of the armature winding or the 
ri ’ u U pe g 
“es battery's internal impedance), 
In order to obtain the analogous values of 8 for the The rectified voltage Uy can also be expressed in 


mode of sinusoidal current at the rectifier inputs, we make terms of the resistance ry: 
use of the curves of current and voltage (Fig. 1b). The ' 
variation in the voltage u at the rectifier inputs can, with Ua = Tara, (11a) 


this, be presented by the following relationship: where, in accordance with (10) 


u=Eh+(Um,—E)sin4, — 
‘eT. +F- (11b) 
O0<Pcne (n=1,3,5,...). (5) 4 
The power applied to the rectifier bridge equals 
Since |u| =| ug| , then the expression for the recti- 





fied voltage will have the same form as in (5), with the p==UI, (12) 
sole difference thatn=1, 2,3,... By substituting (2) and (3) in (12), and taking (11a) 
According to (5), the mean value of the voltage at the into account, we get : 
rectifier output will equal p=s% Z (13) 
bd uri 
oe eat a \ (Um —k) sin d= We introduce the concept of the equivalent load 
-s impedance re, and express the power P in terms of it 
on Bit- SU Bi). (6) P= I, (14) 
where .* U/ 1. 
The effective value of the voltage at the rectifier By substituting (2) in (14), we get 
input is - 
’ tn , Cana P= 14-5. (15) 
U=)V e+“. (7) "i 


By substituting (6) and (7) in (3) we obtain, after 


some transformations, the expression for the coefficient povesien Ser che equivalent impedence 2, in the form 











B , for the mode of sinusoidal current at the rectifier in- 8 
puts: Wea a .. ©? (16) 
ee Um : = ( Tx | : Figure 2 shows the circuit for a single-phase rectifier 
~ E \2 bridge operating on a counter-emf (a) and its equivalent 
E\3 (! a. (8) circuit with active impedance r,(b), Such a transfiguratio® 
lv } ieee. of rectifiers with counter-emf to an equivalent active im 


pedance significantly simplifies the circuit design if a 


By equating (13) and (15), we obtain the definitive ex 
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Fig. 2. 


magnetic amplifier is connected in series with this im- 
pedance, For example, if one ignores the losses in the 
steel of the magnetic amplifier's cores, then the follow- 
ing equation will be valid for the fundamental harmonics 
of the voltages and currents (in order to simplify the nota- 
tion, the subscripts for the fundamental harmonics have 
been suppressed in this and the following equations): 


UL+PQ try =U, 


where U, is the fundamental harmonic of the effective 
value of the inductive component of the voltage across the 
terminals of the magnetic amplifier's working winding 
we, U is the fundamental harmonic of the effective value 
of the source voltage, r_,, is the resistance of working 
winding w~ of the magnetic amplifier. 

By expressing the voltage and current in terms of the 
magnetic induction and magnetic field strength we obtain 
after some uncomplicated transformations, 


(17) 





B2 2 
a Soe f, (18) 
mn, ' i 
where 
Ur J 
Ba = 4.44fw_S’ Bmo = 4.44jw_S’ 
wd w_] w_U 
l(%+ina) 


Here, / and S are the average length and cross-section 
of the magnetic amplifier’s magnetic circuit and I,,3, 
is the maximum current, corresponding to Uy, = 0. 

Equation (18) is the equation of an ellipse with its 
center at the origin and with semiaxes B,,, and Hy By 
joining such an ellipse with the characteristics B,, = 
= f(H, H_), we easily obtain the magnetic amplifier char- 
acteristics we are sceking, A graphicoanalytic method of 
obtaining magnetic amplifier characteristics is explained 
in more detail in [5], 

The magnitude of the resistance re which enters into 
the expression for the magnetic field strength Hy depends 
on the ratio of the coefficients 6 /B.,. For example, for 
a sinusoidal voltage, according to (4a) and (4b), we have: 
for E/ Um = 0, te = 14; for E/ Um 1,1. = 1.41ry. For 
sinusoidal current, according to (9a), (9b), and (9c), we 
have: for E/ Up, = 0, te = 14; for E/ Um = 1s te = 0.9r 4. 

Thus, in determining the characteristics of a mag~ 
netic amplifier by conjoining the magnetic characteristics 
Bm = f(H, H_) with the equation of an ellipse, (18), three 
limiting cases can be found (Fig. 3): E/ Um = 0, f.e., 











Fig. 3. Working trajectories of a magnetic am- 
plifier with a counter-emf, 


there is a purely active load (r,. = ry =th where r ;is the 
load resistance); E/U,, * 1, i.e., the load is a counter- 
emf, Here there can be either r, = 0.9r4 or te = 1.4184. 

In charging storage batteries, principally starter bat- 
teries, there is a significantly large voltage drop from 
their emf, caused by internal impedances, so that E U,,* 
1, In schemes for controlling the rotational speed of 
motor armatures, there can be, as a function of its mag- 
nitude, 0< E/U,, = 1. 

If the voltage drop across a magnetic amplifier is 
greater than the voltage drop across the input (the mag- 
netic amplifier is operating on the upper portion of the 
magnetic characteristics) (Fig. 3), then r, < ri 

The method just presented for obtaining the desired 
characteristic is convenient in those cases when the sup- 
ply voltage and the load parameters are invariant. 

When the load parameters vary, it is convenient to 
give different values to the resistance r, and the current 
I, and to compute the value of H corresponding to them 
from (19), and the value of B,, from the formula 


VRE TP 
Bn = 4. 44jfw_S : 

By conjoining the thus obtained function 8. * f(H) 
with the characteristics Bin = f(H, H_), we can find the 
magnetic amplifier characteristics being sought. 

In many cases, the magnetic amplifier operates with 
H = const, Such modes of amplifier operation occur, for 
example, in charging starter storage batteries, in control- 
ling the rotational velocity of a motor armature with con- 
stant torque, etc, 

In order to obtain the characteristic of a saturated re- 
actor magnetic amplifier, providing a constancy of charg- 
ing current, it is necessary to know the law of variation 
of the voltage across the terminals of the storage bat- 
tery as a function of the charging time for I = const, L.e., 
it is necessary to have the function Ug = f(t) or E = f(t) and 
the internal resistance r, 








(20) 


Magnetic Amplifier Operation with Three- 
Phase Rectifier Bridge 








Three -phase rectifier schemes are widely used in 
technology, chiefly in power rectification stands, How- 
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ever, there are almost no works in the literature which 
consider the questions of designing, or choosing the param- 
eters for, the magnetic amplifiers used with three -phase 
rectifiers, 

In three-phase controlled circuits with magnetic am- 
plifiers, the electromagnetic processes occur in a more 
complex manner than in single-phase circuits. The cur- 
rents and voltages at the input and output of the rectifier 
bridge are different here, even for an active load. 

A three-phase rectifier circuit, just as the single-phase 


circuits considered earlier, can be represented by an equiv- 


alent circuit, For calculating the circuit elements, it is 
convenient to represent a three-phase rectifier circuit by 
a Y-network, With this, the transfiguration must be so ef- 
fected that, in the equivalent Y-network, the same linear 
voltages and currents occur as in the rectifier circuit, 
while the same power as before acts on the load. 

With an active load, the power of one phase of the 
equivalent Y-network, Poe equals 


Ry =U! 


(21a) 


or 
a 
R =l'n, (21b) 
where Tp is the phase resistance of the equivalent Y-net- 
work, 
In a circuit with an active load, the rectified voltage 
ordinarily equals 


Ug = Tara; 


where rq is the active load impedance, 
By setting U = Up in (3), and by substituting (2), (3), 
and (22) in (21a) and (21b), we get 


(22) 


r 


2 
fr 
P, = ]* 32 ; (23b) 
Equating (23a) and (23b), we obtain, finally, 
rp = = Fa (24) 


u 


Expression (24) is the condition for the transfiguration 
of a three-phase rectifier circuit to an equivalent Y-net- 
work with the resistance of each phase equal tor,. This 
expression is analogous to Expression (16) for single ~phase 
rectifier circuits, 

Figure 4 shows the three -phase rectifier bridge cir- 
cuit with an active load (a) and the equivalent Y-network 
corresponding to this circuit (b), An analogous Y-network 
is also equivalent to a three-phase circuit with zero out- 
put, The difference is solely in the magnitude of phase 
resistance r,,, determined by the coefficients 68; and B,,. 

The values of the coefficients 6; and 6 |, for a three - 
phase rectifier bridge circuit operating on an active load 
are given in Table 1, 
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TABLE 1 
Mode of rectifier operation B; BY | 
Sinusoidal voltage on rectifier 
input 1,227 2,12 
Sinusoidal current on rectifier 
input 1,35 2,34 

















For other forms of load (counter-emf and active- 
inductive load), the coefficients 6; and By are given in 
{1, 6, 7). 

Figure 5 shows the circuit of a three-phase rectifier 
with a magnetic amplifier and an active load (a) and the 
equivalent Y-network (b) (for the purposes of simplifica - 
tion, the functional dependence of the magnetic ampli- 
fier*s induction on current is not shown in Fig, 5b). 





























Fig. 5. 


The following equation is valid for the fundamental 
harmonics of voltage and current for each phase of the 
equivalent Y-network [with the same assumptions as for 


(17) }: 


UL +1?» + fa)? = U5. (25) 


In designing three -phase circuits, one frequently has 
to do with linear voltage sources, By taking this into ac- 
count, and by expressing the voltage and current in terms 
of magnetic induction and magnetic field strength, accord- 
ing to (19), after some uncomplicated transformations, we 
obtain 


3b? 
h2 


, 
mo 


3H? 


|. ——— on 
TF 1. (26) 





For the substituting of (19) in (25), it is assumed that 
U is the linear source voltage and that fe = r,. 

As with (18), (26) is the equation of an ellipse which, 
when joined to the characteristics B,, = f (H, H_), allows 
one to obtain easily the desired magnetic amplifier char- 
acteristics, 
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The coefficients 8; and B,, are different for the modes 
of sinusoidal current and sinusoidal voltage on the recti- 
fier input, However, the ratios of the coefficients, B,/ 6 
for each of these modes are almost identical, being equal 
to 0,579 for the mode of sinusoidal voltage, and to 0,577 
for the mode of sinusoidal current, 

Sinusoidal voltage at the rectifier input occurs only 
for Brn = 0. But, By, * 0 always for magnetic amplifier 
operation, Therefore, in computing the magnetic ampli- 
fier characteristics H = f(H_), the ratio Tp = 0.577tg should 
be used, 

If it is required to determine the load characteristics 
H, = f(H.), then it is necessary to multiply the ordinates 
of the curves H = f(H_) by the coefficient 6;. 

For variations of the load parameters or of the line 
voltage U, it is convenient to compute the magnetic in- 
duction B_, by the formula 





ee V U3 — 31? (ry + tma* 
i V3-4.44jw_S 





(27) 


For operation of the rectifier on a counter-emf, the 
equivalent Y-network in accordance with (11a), will have 
a form analogous to that shown in Fig. 4, The only dif- 
ference will be that, in this case, 


Be 


r. rd. (28) 


Pe 

For operation of a three -phase rectifier on a counter- 
emf, the same ratio is retained between the fundamental 
harmonic of the input voltage and the dc component of 
the output voltage as for a purely active impedance, i.e., 
for charging storage batteries with rectified current, the 
coefficients 8 given in Table 1 are valid, Also retained 
is the ratio between the fundamental harmonic of input 
current and the dc component of the output current [1]. 
As for an active load, we consider that Tp = 0.57754, and 
t% is defined by (22), where rg = rm 


Schemes for Devices for the Automatic 
Stabilization of Charging Current 








Charging of many types of storage batteries is car- 
ried out with de which is invariant in magnitude, In 
many case, the invariance in current magnitude is nec- 
essary in the charging process in order that the battery 
charge can be gauged, For example, the charge on acid 
storage batteries is determined by the density of electro- 
lyte and voltage across the terminals, both invariant over 
several hours, for one and the same charging current. The 
charge of alkaline batteries is determined by the magni- 
tude of the capacitance obtained by them (with an invari- 
ant charging current, the charge is here determined actual- 
ly by the batteries’ charging time). 

In charging devices with saturated-core magnetic 
amplifiers, it is not difficult to implement automatic 
stabilization of the charging current (shown in Figs. 6 
and 7), 
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Fig. 7. Fig. 6. 


As a function of the quality of the steel and of the 
operational mode of the magnetic amplifier, the control 
winding Wc, in order to obtain a stabilized charging cur- 
rent, can be connected either to the termiuals of the 
storage battery to be charged, to provide positive voltage 
feedback (Fig. 6), or to.a source of invariable voltage, to 
provide constant bias (Fig. 7). 

To explain the choice of one scheme or anvuther for 
magnetic amplifier control, Fig. 8 shows the characteris- 
tics By, = f(H, H_). 

In the process of charging storage batteries, their 
voltages is changed from Uj, to U% by the amount 


AU, =U4a— Va. (29) 
This change of voltage, in the circuit with positive 

feedback (cf,, Fig. 6), gives rise to an increase in the mag- 
netic amplifier's magnetizing current and, consequently, 
the magnetic field strength due to the current in the posi- 
tive feedback winding is increased by an amount 
AH. = AU,. The quantity AH_ is defined by the formu- 
la 


AH = H..— H. 


(30) 











Fig. 8. Working trajectories of a 
magnetic amplifier used in charg- 
ing devices, 
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where I'_ and H'_ are the magnetizing current of the feed- 
back winding and magnetic field strength corresponding 
to it, both taken at the beginning of the charging of the 
storage battery, and“ and H% are the corresponding 
quantities at the end of the charging. 

Change of the magnetic field strength by the amount 
AHL. gives rise to a decrease in the magnetic induction in 
the magnetic amplifier’s cores by an amount 


ABm = Bu — Bin (31) 
Where B’ andB'',,are the magentic inductions at the 
beginning and end of the charging period respectively. 

In charging storage batteries, their voltage Uy or- 
dinarily changes considerably (in acid storage batteries 
of the starter type, the voltage U%, can reach 135% of its 
initial value U"y)- However, in accordance with (20) or 
(27), to obtain constancy of the charging current, it is 
necessary to decrease insignificantly the magnetic induc - 
tion in the magnetic amplifier's cores, Therefore, the 
magnetic amplifier must operate in the saturated region 
of the characteristic B,, = f(H, H_) so that, with an in- 
crease in voltage Uy, the relative growth in the magnetic 
field strength AH-, induced by the positive feedback cur- 
rent, will be larger than the relative decrease in magnetic 
induction AB,,. In Fig. 8, the working trajectories of the 
magnetic amplifier are shown by the dashed lines between 
points a and b, 

To obtain a charging current which is invariant in 
magnitude, it is necessary to choose such a mode of mag- 
netic amplifier operation in which, for a variation in mag- 
netizing magnetic field strength defined hy (30), the mag- 
netic induction must be changed by the amount AB__, com- 
puted in accordance with (20) or (27) and (31), This mode 
of magnetic amplifier operation, providing automatic stab - 
ilization of the charging current, we express by the cqua - 
tion 

AU, us 


(H. — Hu) neconst = , 1’ 
c 





(32a) 


where (H®_~H'_)}- const 4Fre the limits of variation of the 
magnetizing strength on the family of characteristics B,, = 
= f(H, H_) for variation of the magnetic induction, com- 
puted from (20) or (27), from B4,, to Bin: 

The right member of (32a) is the increment of mag- 
netic field strength as the result of a change in voltage of 
the storage battery in accordance with (29). 

If 


AU 4g Ue 
= ae (32b) 
then the charging current will be decreased in the process 
of charging the battery. In order to provide a constant 
charging current in this case, the magnetic amplifier must 
operate in the regions of low field strength H or high field 
strength H_, i.e., the trajectory must be shifted either to 
the left, or down, or down and to the left simultaneously. 
then the charging current will be increased in the process 


(H— — Hu) n= const > 


448 


» , AU, &% 
(H. — H~) y= const < Ye a ? 
Cc 


(32c) 

of charging the battery. In order to provide here a constant 

charging current, the magnetic amplifier must operate 

in the regions of large field strength H and low field strength 
H= , i.e., the trajectory must be shifted either to the right, 

or up, or up and to the right simultaneously. 

Charging devices in which stabilization of charging 
current was provided by the choice of the proper mode of 
magnetic amplifier operation were considered in [3, 4]. 

To decrease the dimensions and weight of the power 
magnetic amplifiers used in charging devices, it is ordin- 
arily recommended that one use cold-rolled steel such as, 
for example, £320, £330, £370, These brands of steel dif- 
fer from transformer steel, not only by their large values 
of magnetic saturation induction, but also by the steep- 
ness of the slope of the characteristics, particularly in the 
case when toroidal cores prepared from strips are used [8), 

If the characteristics Br = f(H, H.) of the magnetic 
amplifiers used'are almost vertical in their unsaturated 
portions then, to provide constancy of charging current, it 
suffices to provide a constant bias using, for this, small- 
size germanium diodes of type DG~-Ts (Fig. 7). Connecting 
a battery in the bias circuit can permit the elimination of 
the effect of oscillations in the supply voltage U on the 
magnitude of the charging current. Since the magnetic 
induction B,, changes within narrow limits during the bat- 
tery charging process, then this method allows the imple- 
mentation of a completely satisfactory stabilization of 
charging current, 


Example of the Determination of the Basic 
Parameters of a Saturated Reactor Mag- 
netic Amplifier 











We earlier considered methods for determining the 
characteristics of magnetic amplifiers operating on a 
counter emf, We now show that, by means of the relation- 
ships already provided, we can also compute the basic 
parameters of a planned magnetic amplifier. 

We consider the sequence of the determinations of 
the basic parameters of a magnetic amplifier intended for 
the charging of 14 starter storage batteries of type 6S TEN- 
140M by a current of invariable magnitude, equal to 10 
amperes, 

In Table 2 (the upper two rows): there are given the 
charging characteristics of one element of a 6STEN-140M 
storage battery, 50% discharged, for charging by a 10 
amp current, 

The rectified voltage will equal 


U, == mnlg , 


where m is the number of storage batteries to be charged 
(m = 14) and n is the number of elements in one battery 
(n = 6), 

To determine whether or not it is possible to charge 
the given quantity of storage batteries, one must com- 
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| Charging current Iq amperes 10 


_ Voltage on a storage battery 
| element, voits 
| Rectified voltage, volts 
Equivalent phase resistance 
Ip, Ohms 


B= 1 (0) 


to 
_ 
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pute Uy inax* For this, it is necessary to know the maxi- 
mum possible voltage on one element, Ug] max- For 
starter batteries, ordinarily Ug) max = 2.7 v. Then, 
Ug max = 226.8 v. 

According to (3) and the data of Table 1, in the ex- 
ample under consideration, if the voltage drop in the récti- 
fier is neglected, the rectified voltage can equai 
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+e 


V3 





Us= BLU p == (2.12 — 2.34)- = 269 — 297. 


With the given relationships of the possible rectified 
voltage and the maximum possible voltage on the storage 
batteries, we decide that charging will occur with an ap- 
proximately sinusoidal voltage on the rectifier bridge's 
input, But since B,, # 0 with this, then the equivalent 
phase resistance will be defined by the relationship 


o. © 
» = 0.577 jo: 


The results of computing r, are given in Table 2, 

The current I on the constllan input, according to (2), 
will equal 8,15 amp. 

In (27), which relates the given parameters with the 
magnetic amplifier’s parameters, there are still four un- 
known quantities; By,,w~, S and tpyq. The active imped- 
ance of winding w,, is very small and, for a preliminary 
estimate of the parameters, it may be neglected, 


When positive voltaye teedback is used, constancy of 
charging current can be attained only when the magnetic 
amplifier operates on the saturated portion of character- 
istics By, = f(H, H.). Therefore, the magnitude of B,, 
can be given for a well-known brand of steel, By being 
given the magnitude of S as well, we can compute the 
number of turns of winding w.. 

In order that the initially chosen values of w,, and S 
correspond closely to the actual magnitudes, we can use 
the following assumptions, 

For the initial and final moments of storage battery 
charging, the following relationships are valid: 


U? == UG + Pre, U), = 4.44fw SB), 10-8, 
2. yt 4 Dr's s Uy = 4.44fw_ SB", + 10-8, 
Iw 
WW" =W=——- 


l 


The quantities with one prime (*) relate to the be- 
ginning of the charging period of the storage battery, 
while those with two primes relate to the period's end, 
The resistances r,,,, are ignored here, since we are mak- 
ing only a preliminary choice of parameters, 

By transforming these equations, we obtain a func - 
tion for the determination of the volume of steel 





2 _ »2 
vo PV r p—"'p 
4,44fH -10-° V B®, — Be? 


Due to the increase of voltage on the storage batteries 
to provide I = const, the magnetic induction B,,, must be 
decreased, This decrease in induction is ordinarily in- 
significant in the charging of starter storage batteries, By 
taking this into account, and by selecting the magnetic 
induction By, one can compute the value of the square 


root V B®, —B°? 

By being given the magnitude of the magnetic field 
strength H for the selected type of steel, we can approxi - 
mately determine the core volume. 

If the invariance of the charging current is to be pro- 
vided by biasing, then the chosen magnitude of H must 
correspond to the vertical segments of the steel's magnetic 
characteristics, 

The expression given for the volume of steel can 
facilitate to a significant degree the determination of the 
magnetic amplifier parameters if one uses, for example, 
standard three -leg laminas [9]. From the known volume 
V one can choose the type of lamina and, from this, one 
can select / and S, From the known value of S, one can 
approximately determine the number of turns of winding 
w.. by the formula 

IV ip —p 
4.44/S -10-* V B 

From the number of turns found for winding ww, one 
can approximately gauge the resistance r,,,,, for which it 
is necessary to be given the current density and to com- 
pute the mean length of the turns, 

In thus computing approximately the basic param - 
eters of the magnetic amplifier, it is necessary to deter- 
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mine the law of variation of variation of magnetic in- 
duction By = f(t) in accordance with (27) and to con- 
join the function thus obtained for H = const with the char- 
acteristics B., = f(H, H.). The results of the computation 
of Bin = f(t) can be found in Table 2, 

The choice of the required mode of magnetic am- 
plifier operation to provide 1, = const was presented in a 
previous section, 

Limitations on the size of this paper prevent us from 
giving here the entire design of the magnetic amplifier. 


CONCLUSIONS 


The transfiguration suggested of one-phase and three - 
phase rectifier bridges to equivalent loads permits the use 
of the simplest, most graphic and sufficiently accurate 
graphicoanalytic method of designing saturated reactor 
magnetic amplifiers, 

In this paper we have shown the simplest methods of 
stabilizing charging current, which are attained by the 
proper choice of operating mode of the magnetic amplifier. 
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A method is described for determining the mean value of the square of the output quantity of a linear sys- 
tem with n stationary random input processes by means of a simulating device without the use of random func- 
tion generators, This method is a generalization to systems with n inputs of the method described in [2]. 


In solving various problems in the domain of regula - 
tion and automatic control, one frequently has to deal 
with the problem of determining the average value of the 
square of the output quantity of a system with n stationary 
random processes at the system’s inputs, We shall assume 
that the input processes are stationary and ergodic, and 


that their spectral densities (proper and mutual) are rational 


fractional functions of the complex frequency s. 


Posing of the Problem 





We denote by y(t) the output quantity of the linear 
system under consideration. The mean value of the square 
of y(t) is defined by the relationship 


t 

ee ; 

y’ = him Sty *dt. (1) 
0 


A generator of n stationary random processes with 
given proper and mutual spectral densities was described 
in [1]. On the basis of this generator we can construct a 
block schematic for determining the mean value of the 
square of the system's output quantity (Fig. 1). On this 
block schematic, Q;, Q2,... » Q, are the noncorrelated 
generators of white noise, Sh are shaping filters, S,, is the 
given linear system with n inputs, Sq is a squaring device, 
I is an integrator, x,(t), X(t), . . «» x,,(t) are the random 
input processes of system S,, and y(t) is the system’s out- 
put process, 


1,(t) 





Fig. 1, 


For sufficiently large t, one can determine the mean 
value of the square of y(t) from the integrator’s output 
quantity in accordance with (1), With this, it is necessary 
to use n_ generators of noncorrelated white noise, These 
generators must be quite accurate and stable, 

A method was described in [2], applicable to linear 
systems with one stationary random input process, for 
determining the mean value of the square of the output 
quantity by means of a simulating device, eliminating 
the necessity of using white noise generators, This meth- 
od is based on the possibility of generating an input ran- 
dom process with a known spectral density by means of a 
white noise generator and a linear shaping filter Sh (the 
determination of whose transfer function is described in 
(1, 2,]). 

We denote by W(t) the impulsive response of the 
series-connected shaping filter Sh and the given linear 
system S;. Based on known relationships, we obtain 


Loe] 
v= { (IW (IP. (2) 
0 
The block schematic shown on Fig. 2 follows from 
this, In this case, the input of the shaping filter has ap- 
plied to it, not white noise, but a unit pulse function (a 
Dirac function), With increasing t, the integrator's out- 
put quantity tends to the quantity 
The superiority of this method is primarily that the 
white noise generator is excluded from the solving block 
schematic, The linear system is acted upon by a precisely 
defined function of time, which increases the accuracy of 
the determination of the quantity y. Below we give a 
generalization of this method for a linear system S, with 
n stationary random process at the system's inputs, 
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Fig. 2. 
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Solution of the Problem 





The correlation function of the output quantity of 
a linear system with n stationary random processes at the 
inputs is defined [2] by the relationship. 


* by 
meo= Ly Yl) — 


m=] k=10 


— [Wan (41) | Wy (42) Png (F — 42+ 1H) dz,| dts, 
0 


(3) 
where ¢,,,,(T ) is the mutual correlation (cross-correlation) 
function of the mth and kth random processes at the input, 

W,(T ) is the system's impulsive response to the ith input 
on ¢.,,(T) is the correlation function of the system's out- 
put quantity. 

We use (3) to determine the magnitude of y. Since 
the series connection of the shaping filters Sh and the sys- 
tem Ss, is a linear system with n inputs then, in our case, 
W,(T ) will be the impulsive response to the ith input of 
shaping filters Sh at the output of the system S,. The 
function ¢ (7 ) (the correlation functions of the random 
processes at the inputs of shaping filters Sh) are defined by 
the equations 

Pmxk (t) =8(t) for m =k, 
Pmx (t) = 0 for mk. (4) 


By substituting (4) in (3), for T = 0, we obtain the 


quantity y* in the form 
on s ” 
¥* = Fy, (0) = W,(2)]? de. 
w 2} or cores a 
We denote by al 
iti co 
y2 
=\ [W,(*)}? at, 
v= S 1M ™ 


i.e., the mean value of the square of the output quantity 
y_ when the inputs of the shaping filters are acted upon 
by only one white noise generator Q;. 

Then, from (5), we get 


n 
y= 2 #. 
tn (7) 

Consequently, each value of yp {cf., (6)] can be de- 
termined by the method described in [2]. 

Thus, the mean value of the Cone of the output 
quantity of linear system S,, with n stationary random in- 
put processes can be obtained at the output of the block 
schematic shown on Fig, 3, At the inputs of shaping fil- 
ters Sh there act, at the individual moments of time Tj, 
unit pulse functions 5(t-T p- The moments T,(i = 1, 2, 

«+ +» » MN) must be so deomnhtned that, at these moments,the 
entire system is found in a yi gu resting state, The in- 
tegrator I sums the values of y* y step by step ih accordance 
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Fig, 3 





with (7), such that, for t > max { Ty Te ° 
output quantity tends to the value y*, 

It is well known, from the methods of solving problems: 
on simulating devices (analog computers), that the unit - 
pulse functions can be replaced by the introduction of the 
proper initial conditions, 


® Ta} ® its 


CONCLUSIONS 


In the solution of various practical problems, it is 
necessary to determine the mean value of the square of 
the output quality of a linear system S,, at the n inputs 
of which there act n Stationary random processes with 
known spectral densities in the form of rational-fractional 
functions of the complex frequency s. This problem can 
be solved by using the block schematic shown in Fig, 1. 

For a system S, with one stationary random process 
at the input, a method is described in [2] which eliminates 
the use of a random process generator, In the present 
paper, this method was generalized to an arbitrary linear 
system S, with n inputs, acted upon by n stationary ran- 
dom processes with known proper and mutual spectral 
densities, 


The block schematic for the analog computer com- 
putations consists, in addition to the given linear system, 
of shaping filters ( their transfer functions are determined 
by the very simple methods described in [1]) a squaring 
device and an integrator (Fig, 3), At the inputs of the 
shaping filters appear, not white noise, but unit pete, 
functions, With this, the mean value of the square y* is’ 
obtained at the integrator’s output, 

The advantage of this method of computing the 
quantity y“ is, above all, that random function generators 
are eliminated, Since the entire system is acted upon, 
not by random processes, but by accurately defined func- 
tions of time [the pulse functions 5(T )},one can expect 
a greater accuracy in determining the desired quantity 


y. 
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BUREAU OF THE PRESIDIUM OF THE ACADEMY OF SCIENCES OF THE USSR 


Decision of Feb, 12, 1960, Number 134, Moscow 


On the Decision of the Conference on the Theory of Invariance and its Application 





to Automatic Devices of October 20, 1958 (Kiev) 





After hearing the communication of the chairman A. A. Dorodnitsyn, of the commission of the Presidium 
of the AN SSSR on the decision of the Conference on the Theory of Invariance and its Application to Automatic 
Devices, of Oct, 20, 1958 (in Kiev), the Bureau of the Presidium of the Academy of Science of the USSR has de- 


cided; 


1, To endorse the conclusions of the commission of the Presidium of the AN SSSR on the decision of the 
Conference on the Theory of Invariance and its Application to Automatic Devices, of Oct, 20, 1958 (in Kiev) 


the Appendix). 


2. To publish the conclusions of the commission in the journal Avtomatika { Telemekhanika (Automation 


and Remote Control). 


Vice-President of the 
Academy of Science, USSR 
Academician A. V. Topchiev 


Chief Scientific Secretary 
Presidium of the Academy of 
Science, USSR 

Corresponding Member AN SSSR 
E. K. Fedorov 


APPENDIX TO THE DECISION OF THE BUREAU OF THE PRESIDIUM AN SSSR OF 


FEBRUARY 12, 1960, 


Conclusions of the Commission in Connec- 
tion with the Discussion on the Theory of 
Invariance 











The commission, consisting of Academician A. A. 
Dorodnitsyn, Academician AN Ukrainian SSR A. Yu. 
Ishlinskii and Corresponding Member AN SSSR B. N. Petrov, 
designated on October 28, 1958 by the vice president of 
the AN SSSR, Academician A. V. Topchiev, in consider- 
ing the decisions of the Kiev Conference on the Theory of 
Invariance and Its Application to Automatic Devices ( Oct. 
20, 1958), the decisions of the Presidium of the Academy 
of Science of the USSR of April 1, 1941, Section 5, the con- 
Clusions of the commission of the Presidium of the Academy 
of Science of the USSR for discussfon and judging the work 
of Professor G. V. Shchipanov, “Autocontrolled systems 
with several degrees of freedom,” the work of Prof. G. V. 
Shchipanov itself, and certain materials in the discussion 
of the work of G, V. Shchipanov, established the follow - 
ing. 

In the journal Avtomatika i Telemekhanika (Automa- 
tion and Remote Control) No, 1 (1939) there was published 
the paper of Prof. G.V. Shchipanov, “Autocontrolled sys- 
tems with several degrees of freedom”. This paper gave 


Number 134 


rise to wide discussion, the materials of which were pub- 
lished in a number of journals, 

By decision of the Presidium of the Academy of Sci- 
ence of the USSR, taken on March 4, 1940, a commission 
was created to discuss and judge the work of Prof. G, V. 
Shchipanov, “Autocontrolled systems with several degrees 
of freedom ,” this commission consisting of Vice-President 
of the AN SSSR Academician O. Yu. Shmidt, Academician 
S. A. Chaplygin, S. L. Sobolev, N. E. Kochin, and Corre - 
sponding Member AN SSSR N. G. Bruevich, 


The commission's conclusions were discussed at the 
April 1, 1941 session of the Presidium of the AN SSSR, at 
which time the decision was taken to take the commision‘s 
conclusions under consideration, to publish them, and to 
consider the work of the commision as finished, 

The commission's conclusions contained the dissent - 
ing opinions of Academicians V. S. Kulebakin and N, N. 
Luzin which considered _ necessary further research into 
the conclusions of Section 2, in which the estimate of the 
work of G. V. Shchipanov was contained, 

In the following years, a number of works were pub- 
lished which, in some degree or other, treated the questions 
of compensation or invariance (the works of Academician 
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N. N. Luzin, Academician V. S. Kulebakin, A. G, Ivakh- 
nenko , B. N. Petrov, G. M. Ulanov, and others ). 

On October 16-20, 1958, in Kiev, there was held the 
Conference on the Theory of Invariance and Its Applica - 
tion to Automatic Devices, convened by the Engineering 
Section of the AN Ukr SSR, the Kiev City Seminar and the 
Electrical Engineering Institute of the Ukr SSR, in which 
scientific workers from Moscow, Leningrad, Kiev and other 
cities of the Soviet Union participated, 

At the Kiev Conference, a decision (resolution) was 
taken, noting the necessity of developing methods of dis- 
turbance compensation and the principle of invariance, in 
which it was stated, in Section 7: 

"7, In view of the fact that the decision of the Presid- 
ium of the Academy of Science of the USSR on April 1, 
1941 concerning the assertions of the conclusions of the 
special commission on the work of G, V. Shchipanov de- 
layed development of the theory and practice of construct- 
ing automatic systems using compensation of disturbing 
stimuli and various invariance principles and does not 
make possible the proper and necessary progress in this 
important area of automation today, we ask the Presidium 
of the Academy of Science of the USSR to reconsider its 
decision as concerns the erroneous conclusions of the 
aforementioned commission.” 

As a result of the consideration of all the material 
cited above, we can draw the following conclusions: 

1, The conclusions of the aforementioned commission 
to discuss and judge the work of Prof. G. V. Shchipanov, 
taken into consideration and published by the Presidium 
of the Academy of Science of the USSR, must be consider- 
ed correct as concerns the first point and the first paragraph 
of the second point, in which it is stated that: 

"1. The commission considers that questions of con- 
trol theory have undoubtedly tremendous importance in 
the matter of the development of our national economy 
and the strengthening of the defense of our country, and 
must be worked out in the Institute for Automation and 
Remote Control of the Academy of Science of the USSR. 

2. The work of Prof. G. V. Shchipanov under con- 
sideration, "Autocontrolled systems with several degrees 
of freedom,” contains a whole series of incorrect and un- 
founded assertions,” : 

Indeed, G. V. Shchipanov's work did contain incor- 
rect and unfounded assertions, Among these assertions 
were the author's conclusions that, if the compensation 
condition held, then an ideal and universal controller was 
obtained; that the construction of regulators which did not 
satisfy the theorems presented in the paper must be con- 
sidered as unsuitable; that with the compensation con- 
dition met, and for zero initial conditions, a free compo- 
nent appears, etc, 

As is clear from the test and diagrams of G. V. 
Shchipanov's paper, he considered controllers based on 
the principle of deviation without additional direct or 
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indirect stimuli proportional to the disturbance f(t) acting 
on the object of control, and he attempted to use, for the 
synthesis of such controllers, the condition of complete 
compensation while, in this case, it is impossible to use 
the condition of complete compensation. 

2, However, the commission's conclusion (second 
paragraph of the second point) that the basic error in Prof, 
G, V. Shchipanov's paper is the "compensation condition®, 
which leads mathematically to absurd conclusions and 
which does not correspond at all to actually realizable 
controllers, contains itself an erroneous conclusion, in 
that it is far too general a statement. 

The basic error of G. V. Shchipanov consisted, not in 
the formulation of this condition, as was asserted in the 
commission's conclusions (on the contrary, this was his 
fundamental result!), but in the application of the com-~ 
pensation condition to the design of a class of control sys- 
tems which does not permit its use, 

As was shown by subsequent investigations, including 
those in Proceedings of the Kiev Conference on Invariance 
Theory, “the compensation condition” or the "invariance 
condition” can be used for designing a number of types of 
automatic control systems, for example, combined systems, 
with stimulus by deviation and by disturbance, intercon- 
nected automatic control systems for several quantities, 
and also several types of dynamic systems which do not 
come under the heading of automatic control systems, for 
example, bridge systems, complex electric networks, gyro- 
scopic devices, etc. 

Thus, the "compensation condition", or the “invarti- 
ance condition" formulated by Prof. G. V. Shchipanov was 
a new design relationship which can be successfully used 
for designing a definite class of dynamic systems, 

3, The aforementioned inaccuracy in the formula- 
tion of the commission's conclusions, and the absence of 
constructive criticisms in them, led, in some cases, to in- 
correct judgments on works in automatic control, even 
when the compensation condition was used correctly in 
them, which could not render possible development in * 
this direction, Therefore, the commission considers it 
advantageous to publish, in one of the technical journals 
an expanded paper devoted to the prospective applications 
of the invariance principle, and which would cite the 
cases when it is advantageous to use this principle and 
when its use is inadmissible, 

At the same time, the commission considers it nec- 
essary to note the imperfect editing of the decision of the 
Kiev Conference on the Theory of Invariance and its Ap- 
plication to Automatic Devices, from which one might 
gather the impression that the works of G. V. Shchipanov 
did not contain errors, 

Commission Members 
A. A. Dorodnytsin, 
A. Yu, Ishlinskii, 
B. N, Petrov, 





